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The change of hadron properties in dense matter based on various theoretical approaches are 
reviewed. Incorporating these medium effects in the relativistic transport model, which treats 
consistently the change of hadron masses and energies in dense matter via the scalar and vector 
fields, heavy-ion collisions at energies available from SIS/GSI, AGS/BNL, and SPS/CERN are 
studied. This model is seen to provide satisfactory explanations for the observed enhancement of 
kaon, antikaon, and antiproton yields as well as soft pions in the transverse direction from the SIS 
experiments. In the AGS heavy- ion experiments, it can account for the enhanced K + /ir + ratio, the 
difference in the slope parameters of the K + and K~ transverse kinetic energy spectra, and the 
lower apparent temperature of antiprotons than that of protons. This model also provides possible 
explanations for the observed enhancement of low-mass dileptons, phi mesons, and antilambdas in 
heavy-ion collisions at SPS energies. Furthermore, the change of hadron properties in hot dense 
matter leads to new signatures of the quark-gluon plasma to hadronic matter transition in future 
ultrarelativistic heavy-ion collisions at RHIC/BNL. 
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I. INTRODUCTION 



An important feature of quantum chromodynamics (QCD) is its approximate SU(3)l x SU(3)r chiral symmetry 
as a result of small quark masses [[lj. This symmetry is, however, spontaneously broken, leading to a finite quark 
condensate in vacuum, (qq)o ~ (uu)o ~ (dd}o, as given by the Gell— Mann-Oaks-Renner relation, 

m 2 jl = -2m q (qq) . (1) 

In the above, m, ~ 138 MeV and /„• » a^93 MeV are the pion mass and decay constant, respectively; and m q = 
(m u +rrid)/2 ~ 5.5 MeV is the average up and down quark masses. The quark condensate in vacuum thus has a value 
(qq)o ~ -(245 MeV) 3 . 

As the density and/or temperature of a hadronic system increase, this spontaneously broken symmetry is expected 
to be partially restored, so the quark condensate would decrease with increasing density and/or temperature. At zero 
baryon chemical potential, the decrease of quark condensate with increasing temperature has been observed in lattice 
QCD simulations || as well as in calculations based on the chiral perturbation theory || and the interacting pion 
gas model 0. At finite baryon density, model- independent studies using the Feynman-Hellmann theorem J^|^] have 
shown that the ratio of the quark condensate in medium to its value in vacuum is given by 

(qq)o ftm* 

where pn is the nuclear density and S^at ~ 45 MeV is the ttN sigma term. At normal nuclear matter density po w 0.16 
fm~ 3 , the condensate is seen to decrease already by s» 1/3. Higher-order contributions to the quark condensate in 
nuclear matter due to nucleon-nucleon (NN) interactions have been studied using various models |j]-[nj. It is found 
that at normal nuclear matter density they change the leading-order result by only about 5% [6 11 1. The decrease 



of quark condensate in medium may lead to reduced hadron masses as shown in QCD sum-rule studies |12| , |13| . The 
study of chiral symmetry breaking and partial restoration in hot dense matter is a topic of great current interest in 
nuclear physics [jl4| ^6|. 

In high energy heavy-ion collisions, theoretical simulations based on transport models have shown that a hot dense 
hadronic matter is formed in the initial stage of the collisions. This provides thus the possibility of studying the 
partial restoration of chiral symmetry through the changes of hadron properties in medium |l^-|l^| . Unfortunately, 
the dynamics of heavy-ion collisions is very complex, involving a violent initial compression, which is then followed 
by a relatively slow expansion and finally reaches the freeze out when particle interactions become unimportant. The 
entire reaction typically lasts for about a few tens fm/c. The interesting physics of hadron in-medium properties and 
chiral symmetry restoration can only be studied for a few fm/c during the early part of the expansion stage when both 
the density and temperature of the hadronic matter are high. This stage of the collision can in principle be probed 
by detecting the emitted electromagnetic radiation such as the real and virtual (dilepton) photons. However, both 
are not easy to measure due to their small rates. Furthermore, they can also be produced from initial hard collisions 
and final hadron decays, and this makes it difficult to extract the signals from the hot dense matter. What are 
usually measured in heavy-ion experiments are instead the momentum distributions of hadrons, such as the nucleon, 
nuclear clusters, pion, kaon, etc., which are mostly ejected from the colliding system at freeze out. To infer what have 
happened in the initial hot dense matter from the final hadron phase space distributions requires thus a model that 



2 



can describe the whole collision process. Indeed, various transport models have been developed during the past ten 
years for this purpose. We shall review briefly the history of the development of transport models. We note that a 
number of good review articles are available on both relativistic and nonrelativistic transport models for heavy-ion 
collisions at various energies p(i|-p7j|. 

Since the reduced hadron masses in medium can be consistently included in the relativistic transport model |2§|-|3"c|] 
we shall concentrate mostly on the results obtained from this model. In the relativistic transport model, which is 
derived from the Walecka model the effective hadron mass is connected to the scalar field while its energy 

is shifted by the vector potential. Because of dropping hadron masses due to the attractive scalar field, particle 
production is enhanced as a result of reduced threshold and increased phase space. Furthermore, the dense nuclear 
matter also provides a strong mean-field vector potential, which then affects the momentum distributions of hadrons. 
The relation of the Walecka model to the chiral effective field theory has recently been explored in Refs. [^3|-^5|. 
The attractive scalar and repulsive vector potentials for a nucleon in nuclear matter have also been obtained in the 
Dirac-Brueckner-Hartree-Fock (DBHF) approach |||-[4l| as well as the QCD sum-rule studies [H[|2[Q. 

For heavy-ion collisions at energies available from the SIS at GSI, the relativistic transport model gives a satisfactory 
explanation for the observed enhancement of kaon p5|-^8[ , antikaon (4^j5(| , and antiproton yields |^l],[52| as well as low 
energy pions in the transverse direction p3|-|55[| . Experiments on kaon flow have also provided important information 
on the kaon potential in dense matter |5(],^7j . Furthermore, the study of dilepton production from these collisions is 
expected to allow us to probe directly the in- medium properties of vector mesons in dense matter |58|]59(| . 

The relativistic transport model has also been used to describe the expansion stage of heavy-ion collisions at energies 
available from the AGS at BNL. It can account for the enhanced K + /n + ratio | |60| |, the difference between the slope 
parameters of the K + and K~ transverse mass spectra JtUj ], and the lower apparent temperature of antiprotons than 
that of protons |32| as observed in experiments |63|,|64| . 

For the expansion stage of heavy-ion collisions at energies available from the SPS at CERN, it has been found 
that the relativistic transport model can explain quantitatively the observed enhancement of low-mass dileptons if 
one includes the decrease of vector meson masses in hot dense matter [j35l-|69|l . In a simplified hydrochemical model, 
medium effects also lead to enhanced production of phi mesons [[70 71 1 and antilambdas |^,^3) in these collisions. 

The success of theoretical studies including the medium effects in explaining a large body of experimental data 
from high energy heavy-ion collisions thus provides the possible evidence for the partial restoration of chiral sym- 
metry in hot dense matter formed in high energy heavy- ion collisions. The purpose of this article is to review the 
theoretical predictions for hadron properties in dense matter, the consistent incorporation of these medium effects in 
the relativistic transport model, and the observable consequences in heavy-ion collisions from SIS to SPS energies. 

In Section II, hadron properties in dense matter are discussed using various theoretical approaches, such as the 
Walecka model, the DBHF approach, the QCD sum-rule approach, and the effective hadronic model. The relativistic 
transport model based on the nonlinear a-uj model is then described in Section III, together with a brief review of the 
development of transport models. In Section IV, the results obtained from the relativistic transport model for particle 
production and flow in heavy-ion collisions at SIS energies, the strangeness enhancement in heavy-ion collisions at 
AGS and SPS energies, and dilepton production in heavy-ion collisions are then presented and compared with the 
experimental data. The relevance of these experimental observables as the probes of hadron properties in dense matter 
is emphasized. A summary and outlook is then given in Section V. 
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II. HADRONS IN DENSE MATTER 



Theoretically, lattice QCD simulations should provide us with the most reliable information about the temperature 
and/or density dependence of hadron properties. Currently, these simulations can be carried out only for finite 
temperature systems with zero baryon chemical potential. Furthermore, the quenched approximation is usually 
introduced for quarks, so their dynamical role is only crudely included. From the hadronic correlation functions at 
large spatial separations the hadron screening masses have been obtained in the lattice QCD f74|-|76||. In a recent 
paper, Boyd et al |f7^| have shown that, up to about 0.92T C , where T c is the critical temperature for the quark-gluon 
to hadronic matter phase transition, there is no significant change of the rho-meson screening mass. They have also 
found that the quark condensate does not change very much until one is extremely close to the critical temperature. 
This is in contrast with the result from chiral perturbation theory which shows a more appreciable decrease of the 
quark condensate with temperature Q. Therefore, the study of hadron properties in dense nuclear matter at present 
time has to rely on theoretical models. In this section, we shall review some of these approaches. In particular, the 
properties of baryons (mainly nucleon), pseudoscalar mesons (mainly pion and kaon), and vector mesons (mainly rho, 
omega, and phi mesons) will be discussed. 

A. baryons 

A simple model for nuclear matter at high density is the Walecka model, also known as the quantum hadrodynamics 
(QHD-I) or linear a-u model, in which nucleons interact with each other through the exchange of a scalar sigma and a 
vector omega meson pl| , p2| . The sigma meson gives rise to an intermediate-range attraction while the omega exchange 
leads to a short-range repulsion, both are the essential properties of nuclear force. In this phenomenological model, 
the nucleon mass is reduced by the attractive scalar interaction, while its energy is shifted by the repulsive vector 
interaction. The parameters in this model are adjusted to fit the nuclear matter properties at saturation density using 
the mean-field approximation. 

This model gives, however, too small a nucleon effective mass and too large a compression modulus at saturation 
density compared to what are known empirically. To correct for these deficiencies, self interactions of cubic and 
quartic forms have been introduced for the scalar field |7^,f78| . The Lagrangian for this nonlinear a-u> model is given 
by 

C = N[^(id» - g^ 1 ) - (m N - g a N)]N + ^(d^d^a - m 2 ^ 2 ) - ^ba 3 - ^ca 4 
- 1(0^ - d v Lj^(dW - d v ^) + ^mX< (3) 

where N is the nucleon field with mass to at, while a and w M are the fields for the scalar and vector mesons with 
mass m a and m w , respectively. The parameters b and c determine the strength of scalar field self- interaction. With 
b = c = the nonlinear a-uj model reduces to the original Walecka model. The coupling constants between a nucleon 
and the scalar and vector fields are denoted by g a and g u , respectively. In a static nuclear matter, meson fields in 
the mean-field approximation become constants and are independent of the spatial coordinates. Furthermore, the 
vector field has a nonvanishing value only in its time-like component due to translational invariance in infinite nuclear 
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matter. The equation of motion for a nucleon with momentum p in a nuclear matter of density pjv is then given by 
the Dirac equation 

(a-p + Pm* N )N = e*(p)iV, (4) 
where the nucleon effective mass is related to the scalar potential S5 = —g cr (a), i.e., 

m* N = m N + Eg, (5) 

2 

and its energy is shifted by the vector potential £y = ^-p_/v, i.e., 

£ * = (m^+p 2 ) 1/2 + £°v. (6) 

While the vector potential is simply proportional to the nuclear density pn, the scalar potential depends on the 
scalar density ps via the following self-consistent equation, obtained by minimizing the energy, 

ml(<7)+b(a) 2 + c(<j) 3 =g aPs . (7) 



Table I Parameters of the nonlinear a-tu model corresponding to the soft and stiff equations of state and the nuclear 
matter properties at the saturation density po = 0-16 fm -3 . 

C = S^mjv C = g^rn^ g _ (J = _^ S ( M mfc R ( M y) 
a ma " m, 9%m N g% A v ' m N V '_ 

soft 13.95 8.498 0.0199 -0.00296 -15.96 0.83 200.0 

stiff 15.94 12.92 8.0 x 10~ 4 2.26 x 10" 3 -15.96 0.68 380.0 



The four parameters in the nonlinear a-u model are determined by the nuclear matter saturation density (po=0.16 
fm -3 ), binding energy (B=15.96 MeV), compression modulus (K), and the nucleon effective mass (m*)- In this review, 
only two sets of parameters, corresponding to the soft and stiff equations of state (EOS), are considered. At po, the 
soft EOS corresponds to m* N /m N = 0.83 and K=200 MeV, while the stiff EOS has m* N /m N = 0.68 and K=380 MeV. 
The corresponding parameters are listed in Table 1. The two nuclear equations of state, together with the nucleon 
effective mass, are shown in Fig. 1 by solid curves. Both the soft and the stiff equation of state in this relativistic 



formulation are very close to that from the Skyrme interaction 24 shown in the right panel of Fig. 1 by dotted 
curves. 

It is worthy to mention that in the literature the expression 'effective mass' has been used to denote different 
quantities. A nice clarification of the relationship between these various quantities can be found in Ref. fnl. In 
nonrelativistic models, the nucleon effective mass is related to the nonlocality, or momentum dependence, of the 
single-particle potential. This quantity should be identified with the Lorcntz mass ]79[ | defined in terms of the 
Schrodinger-equivalent potential in relativistic models, rather than the Dirac mass defined in terms of the scalar 
potential as in Eq. (||). 

Extensions of the Walecka-type model from zero to finite temperature have been studied in Refs. |80|-|82]]. It is 
found that at finite density the nucleon effective mass first increases slightly at low temperature and then decreases 
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rapidly at high temperature. A similar temperature dependence of the constituent quark mass has been found in 
studies based on the Nambu— Jona-Lasinio model Ml. 




(ASH) ^TJU 



Fig. 1 Left panel: The nuclcon effective mass in medium. Solid and dashed curves are from the nonlinear a-uj 
model and DBHF theory, respectively. Right panel: Nuclear equation of state. Dotted curves are from the Skyrme 
interaction. 



A more microscopic approach to nucleon properties in nuclear matter is provided by the DBHF calculation |37 
using, e.g., the Bonn one-boson-exchange potential |36| , |37| |. As in conventional Brueckner approach, the basic quantity 
in DBHF approach is the G-matrix which satisfies the in-medium Thompson equation |37],|{],[4l| , 



G(q',q;P) 



+ 



V(q',q) 

d 3 k , m 



Q(k,P) 



(2tt) 3 



^"P/2+k -^P/2+q 2-E p / 2+k 



G(k,q;P), 



(8) 



with E£ — (m*^ -fk 2 ) 1 / 2 . The momenta P and q denote, respectively, the center-of-mass and the relative momentum 
of two interacting nucleons. The kernel of this integral equation is given by the sum of the one-meson-exchange 
amplitudes from six nonstrange bosons (n, er, p, uj, rj,S). A form factor of monopole type is introduced at each meson- 
nucleon vertex to take into account the short-range dynamics due to quarks and gluons. The G-matrix becomes 
density-dependent through the Pauli projection operator Q and the scalar field S5, which appears in m* N as shown 
by Eq. (|). 

The Dirac equation and the in-medium Thompson equation are then solved self-consistently to determine the G- 
matrix, from which the nuclear matter properties can be derived. For example, the nuclear equation of state, or the 
energy-per-nucleon, is given by 



£ f x ^ rriNrn^ 



1: 



94 ^ 



'A' 



2 A ^ E*E* 
i,j<k F 1 1 



{ij\G\ij - ji) - m N , 



(9) 



where the nucleon bare mass has been subtracted out. 



G 



The G-matrix also allows one to determine the nucleon scalar and vector potentials in nuclear matter, i.e., 

^ s + ^ v =ReJ2 ^(ij\G\ij-ji)- (10) 

1 j<k F 1 3 

The nuclear equation of state and the nucleon effective mass obtained from the DBHF calculation with the Bonn A 
potential are shown in Figs. 1 by dashed curves. 

Nucleon properties in dense matter have also been studied in approaches based on the QCD. In the QCD sum- 
rule approach ||[l^,[l2-44 8^j8^], hadron properties are related to a number of condensates, which are the vacuum 
expectation values of quark and gluon fields and describe the nonperturbative aspects of the QCD vacuum. The 
QCD sum rules have been quite successful in describing the hadron properties in free space, e.g, the nucleon mass is 
quantitatively reproduced in this approach J86|-p^| . The extension of QCD sum rules to hadrons in nuclear matter is 
achieved by using in-medium condensates. For a nucleon, it has been shown that the change of scalar quark condensate 
in medium leads to an attractive scalar potential which reduces its mass, while the change of vector quark condensate 
leads to a repulsive vector potential which shifts its energy [ fi^|]4^ -|44f| . The nucleon scalar and vector self-energies in 
this study are given, respectively, by 

~ 871-2 II- \ 871-2 ^ N 

Es W "]Sf - {qq)o) " ~M%m^Tm^ 

Mir 2 , , . 32tt 2 . , 

where the Borel mass Mb is an arbitrary parameter. With Mb ~ mjv, and m u + rrid w 11 MeV, these self-energies 
have values which are similar to those determined from both the Walecka model and the DBHF approach. 

Experimental evidences for these strong scalar and vector potentials have been inferred from the proton-nucleus 
scattering at intermediate energies J)0| via the Dirac phenomenology |)1 92 1 in which the Dirac equation with scalar 
and vector potentials is solved. The proton-nucleus scattering has also been successfully described using realistic NN 
interactions in the impulsive approximation |9^j9^| or the local density approximation [p5f . 

In heavy-ion collisions at energies above the pion production threshold, baryon resonances can be produced. Trans- 
port model simulations have shown that in central heavy-ion collisions at incident energies around 1 GeV/nucleon, 
the baryon resonances to nucleon ratio can be as large as 0.2 |45|j . This ratio approaches 1 when the incident energy 
is increased to about 10 GeV/nucleon at AGS pq , ^7| . The mean-field potentials for resonances are thus needed in 
transport models, and up to now have been assumed to be the same as those for a nucleon. This assumption for 
a delta particle seems to be supported by analyses of the pion-nucleus scattering data [^8[ . However, recent QCD 
sum-rule studies [ p9[ have shown that the delta vector potential is considerably weaker than that for a nucleon, while 
its scalar potential is somewhat stronger. This would lead to a larger net attraction for a delta resonance than for a 
nucleon in nuclear medium. 

Also of interest is the potentials for hyperons in nuclear matter. As for nucleon, most of the empirical information 
about hyperons in nuclear matter have been extracted from the structure of hypernuclei. The potential for the lambda 



particle is relatively well determined to be about -30 MeV [ 100 1 . The separation of this nonrelativistic potential into 



the scalar and vector potentials in the relativistic approach is, however, not without ambiguities [ 101 1 . 



The lambda potential in nuclear matter has also been studied using the DBHF approach [102]. The result ranges 
from -25 to -40 MeV, depending on the input boson-exchange models for the AN interaction. There have also been 
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various attempts to generalize the Walecka-type model from SU(2) to SU(3) to include the hyperon degrees of freedom 



[ 103 - 106 1 . In the naive SU(3) quark model, the hyperon potential is about 2/3 of the nucleon potential, as there 
are only two light quarks in a hyperon instead of three light quarks in a nucleon. Recently, hyperon properties in 
nuclear matter have also been studied using the QCD sum-rule approach. It has been found that both lambda scalar 



and vector potentials are significantly weaker than the prediction from the naive quark model [107|, while those of 



the sigma hyperon are somewhat stronger and are close to the ones for the nucleon |108|. The accuracy of these 
findings are, however, limited by uncertainties in the nucleon strangeness content and certain in-medium four-quark 
condensates. 

B. pseudoscalar mesons 

Pseudoscalar mesons play a special role in nuclear physics. On the one hand, they are the Goldstone bosons of 
spontaneously broken chiral symmetry and are thus closely related to chiral symmetry restoration in nuclear medium. 



Furthermore, their small masses make the chiral perturbation theory applicable for low energy phenomena [ 109 - 1 12 1 . 
On the other hand, in heavy-ion collisions, pions are the most copiously produced particles; e.g., at SPS energies, 
the pion/nucleon ratio approaches 10 in midrapidity region. Knowledge on the properties of pseudoscalar mesons in 
medium is thus very important for understanding both chiral symmetry restoration and the dynamics of high energy 
heavy-ion collisions. 



1. pion 



Pion properties in dense matter have been a subject of extensive studies over many years [ 1 13 1 . The strong p-wave 



pion-nucleon interaction through the delta resonance is known to modify appreciably the pion dispersion relation at 



finite density [113 1 15 1 . In the delta-hole model, the pion dispersion relation in nuclear medium can be written as 

io(k,p B ) =ml + k 2 +n(c;,k), (12) 



where the pion self-energy is given by 



n(w, k) = — — — , (13) 

1 - 5 X(w, k) 



with g' « 0.6 being the Migdal parameter due to short-range correlations. The pion susceptibility \ is given by 

x(c,k)«|(^) 2 ^^exp(-2kV6 2 )p., (14) 

9 V 77V / UJ Z — UJ R \ / 

where f^NA ~ 2 is the pion-nucleon-delta coupling constant, b « 7m w is the range of the form factor, and ujr ~ 
j£-+m A -m N . 

The pion dispersion relation obtained from this model is shown in Fig. 2. The pion branch in the lower part of 
the figure is seen to become softened, while the delta-hole branch in the upper part of the figure is stiffened. These 
results are, however, modified by the imaginary part of the pion self-energy due to the delta decay width and pion 
absorption [116- 12(j |. As a result, the strength of the delta- hole branch is significantly reduced, and the peak in the 



spectral function is less shifted in energy than the pion branch obtained without the imaginary pion self-energy. 
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Pion properties in dense matter have also been studied using the chiral perturbation theory |121 122]. In effective 
chiral Lagrangians a pion has a negligible vector interaction with a nucleon. The s-wave irN scalar interaction, 
proportional to the ttN sigma term S^at, would lead to a reduction of the pion mass |l2l| . This effect is, however, 
almost completely canceled by the range term from the energy dependence of the ttN scattering amplitude, so the net 



mean-field effects are quite small 123 , and the pion mass is found to increase only slightly with increasing density. 



This is consistent with the results determined from the pionic atom data [ 1 13 1 . 

The Nambu— Jona-Lasinio model has also been used to study pion (or in general Goldstone boson) properties at 



both finite density and temperature [124,125]. In this model, pseudoscalar mesons are treated as quark-antiquark 
collective modes. Because of the finite scalar mean squared radius of a constituent quark, the s-wave nN attraction 



is found to be screened, leading to a slow increase of the pion mass with increasing density and/or temperature [ 125 1 . 
On the other hand, the pion mass is found to decrease slightly with temperature in an effective Lagrangian that 



includes its interactions in a pion gas via the rho and a\ mesons [126 
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Fig. 2 The pion dispersion relation from the delta-hole model. (From Ref. [115]) 



2. kaon 



The study of kaon properties in dense matter has recently attracted a lot of attention as it is relevant for un- 
derstanding not only the strangeness enhancement in heavy-ion collisions, which has been suggested as a possible 



signature of quark-gluon plasma formation in these reactions [127], but also the kaon condensation in neutron stars 



[128 1 and the possible formation of mini-black holes in galaxies [12£ 
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The properties of kaon in dense matter were first studied in chiral Lagrangian by Kaplan and Nelson |130| , and 
further pursued by many others [ 1H 1121 , 131 138 j. The SU(3)ixSU(3)fl nonlinear chiral Lagrangian used by Kaplan 
and Nelson is written as [130,134 136 

C = i/ 2 Tr^S^E+ + i/ 2 A[TrAf 9 (£ - 1) + h.c.] + TvB(t^d^ - m B )B 
+ iTrBj^V^B] + DTrB^~/ 5 {A^ B} + FTtBj^j 5 [A^ , B] 
+ aiTiB(£,M q i + h.c.)B + a 2 TrB.B(£M 9 £ + h.c.) 

+ a 3 [TrAf 9 S + h.c.]TrBB, (15) 

In the above, B is the baryon octet with a degenerate mass tub, and 

S = exp(2wr//), and £ = VS = exp(«r//), (16) 

with 7r being the pseudoscalar meson octet. The pseudoscalar meson decay constants are equal in the SU(3)y limit 
and are denoted by / = f n ~ 93 MeV. The meson vector and axial vector currents are defined as 

V, = \{ed^ + £cU + ) and A^ = - £d^+), (17) 

respectively. The current quark mass matrix is given by M q = diag{m 9 , m q , m s }, if we neglect the small difference 
between the up and down quark masses. 

Terms involving the axial vector current can be ignored as they have no effects on the kaon mass. Expanding E to 
order of l// 2 and keeping explicitly only the kaon field, the first two terms in Eq. (Fl™ can be written as 



& i Kd li K-A(m q + m s )KK + ---, (18) 

where 

(K+\ 

K = I Q j and K = (K~ K°), (19) 

and the ellipsis denotes terms containing other mesons. 

Keeping explicitly only the nucleon and kaon, the third and fourth terms in Eq. (|15|) become 



3? 

N{i^d„ - m B )N - — N~/°NK d t K + ---, (20) 

where 



N = [ P ) and N = (p n), (21) 



n 



and the ellipsis denotes terms involving other baryons and mesons. 

The last three terms in Eq. (|l^) can be similarly worked out, and the results are 



Tr£(£M g £ + h.c.)B = 2m q NN - — -{m q + m s )KK + 



TrBB(£M q £ + h.c.) = 2m s NN -—(m q + m s )KK + ■ 



NN . 

W 2 

NN 

T 2 



2NN 

[TrM q T, + h.c.]TrBB = 2(2m q + m s )NN T^-{m q + m s )KK H . (22) 
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Combining above expressions, one arrives at the following Lagrangian, 



C = N(i^dn - m B )N + d^Kd^K - A(m 9 + m s )KK 

- ^Nj°NK d t K+ [2m qai + 2m s a 2 + 2(2m q + m s )a 3 ]NN 
NNKK 

2^2— (m q +m s )(a 1 +2a 2 +4a 3 ) + ---. (23) 

From this equation, one can identify the kaon mass by 

m 2 K = A(m q + m s ) (24) 

and the nucleon mass by 

mpf = tub — 2\a\m q + ci2m s + a 3 (2ra q + m s )]. (25) 
Also, the KN sigma term can be expressed as 



^kn = ^{m q + m s )(N\uu + ss\N) 
1 drriM dm,N 



1 

2< 
1 

2 

1 

~2 



2 dm q dm s 



(m q + m s )(ai + 2a 2 + 4a 3 ). (26) 



The second line in the above follows from explicit chiral symmetry breaking in the QCD Lagrangian, and the last 
step is obtained using Eq. (]25|). Then Eq. ( |2^ ) can be rewritten as 

T 2 



C = Niyfdp - m N )N + d^Kd^K - (m\ j£-NN)KK 



^-N^NK d t K+---. (27) 



The above Lagrangian does not describe properly the nuclear matter properties. In Ref. [ 134 1 , an attempt has been 
made to combine this SU(3)lxSU(3).r chiral Lagrangian with the Walecka model |n],[32| by allowing the nucleon to 
also couple with the phenomenological scalar (a) and vector (oj) fields so that the saturation properties of nuclear 
matter can be obtained in mean-field level. Possible justifications of this approach have been further studied in Ref. 



[139 1 based on the chiral perturbation theory. Following this approach, one modifies Eq. ( |27| ) to 

C = Nfr'fdn -m N + g a cr)N - g u N'fNu> lt + C {a, u M ) 

+ d^Kd^K - (m 2 K - —j^-NN)KK - ^N^NK d t K+---. (28) 

In the above, Co contains the kinetic and mass terms of the a and u> fields as well as possible cubic and quartic a 
self-interaction terms. 

From the Euler-Lagrange equation and using the mean-field approximation for the nucleon field, i.e., replacing 
(N"y N) by the nuclear density p^ and (NN) by the scalar density ps, one obtains the following Klein-Gordon 
equation for a kaon in nuclear medium, 



+ ^ PN d t + {ml - ^jfps) K = 0. ( 2!) ! 
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The kaon dispersion relation in nuclear matter is then given by 

^ 2 (k, Pb) 



, 2 ^KN , 3 U! 

k - ~J2~ PS + \Ji pN ' 



(30) 



where k is the three-momentum of the kaon. The third term in the above equation is from the attractive scalar 
interaction due to explicit chiral symmetry breaking and depends on the kaon-nucleon sigma term T,kn- With a 
strangeness content of the nucleon, defined by y — 2(N\ss\N) / (N\uu + dd\N) w 0.1 — 0.2 as normally used, its value 
is 370 < Y>kn < 405 MeV if one takes the strange quark to the light quark mass ratio to be m s /m as 29. On the 
other hand, recent lattice gauge calculations [140 141 show that y « 0.33 which would give T,kn ~ 450 MeV. The 
last term in Eq. (|3^) is from the repulsive vector interaction and is proportional to the nuclear density pn- For an 
antikaon this term becomes attractive due to G parity. 

From the dispersion relation, the kaon energy in medium can be obtained, i.e., 

21 1/2 



aj(k,p N ) 



m 2 K 



J KN 



P 



3pN 

8f 



3 Pn 
8 f 2 ' 



(31) 



Using the KFSR relation [m p = 2^/2 fg p ) and the SU(3) relation (g u — 3g p ) 142 1, the kaon vector potential in the last 
term can be written as (1/3)(<7 w /to w ) 2 /5a/- which is just 1/3 of the nucleon vector potential. This can be understood in 
the constituent quark model as the kaon contains only one light quark as compared to three light quarks in a nucleon. 
Since an antikaon has one light antiquark, its vector potential becomes attractive, leading to a reduction of its energy 
in medium. 



Corrections to the mean-field results of chiral Lagrangian have been studied in Refs. [11C-112]. The term which 
contributes at the same order in chiral perturbation theory as the Kaplan-Nelson term is the energy-dependent range 
term, which can be included by multiplying the Kaplan-Nelson term by a factor 1 — 0.37uj 2 /m 2 K |p3. In Ref. KfJ], a 
smaller value of Tjkn ~ 350 MeV has been used in evaluating the kaon scalar potential so that the effect of range term 
is effectively included. Also, the Nambu— Jona-Lasinio model has been used to estimate these corrections, and they 
are found to cancel off almost completely the attractive scalar potential (TJJ. One notes, however, that the former is 
based on a systematic expansion while the latter is model-dependent. 

The in-medium kaon mass is normally given by the kaon energy at zero three-momentum. With Yjkn ~ 350 MeV, 
the kaon mass increases by about 10 MeV at normal nuclear matter density as shown in Fig. 3. The antikaon mass 
is also shown and is seen to decrease with density as its vector potential is also attractive. 

Because of baryon number conservation, vector potentials in the initial and final states of a hadron-hadron reaction 
are similar and can thus be neglected in estimating the threshold for particle production in medium. It is then 
sometimes useful to consider the kaon effective mass by including only effects due to the scalar field, i.e., 

1/2 



m K 



1 



/H ps 



(32) 



for both kaon and antikaon. In Ref. |136| , the quadratic vector term, ff > m Eq. (|T]) was discarded based on 
the expectation that it would be cancelled by higher order terms not included in the calculation. Then, the kaon 
"effective" mass should involve only the scalar field, while the vector mean field shifts the kaon chemical potential or 
energy, as in Walecka mean-field theory. 



12 



From Eq. ( |3l| ) one can define the kaon potential as the difference between its energies in the medium and in free 
space |5r]Jl43|1 . In the mean-field approximation this potential at normal nuclear matter density is about 10 MeV for 
a kaon at rest. This is smaller than what is expected from the impulse approximation using the kaon-nucleon (KN) 
scattering length in free space which gives about 30 MeV. The inclusion of higher-order terms [ 1 1 0| 1 1 2 1 and possible 



scaling of the pion decay constant / in medium j34 , 144 , L45 would bring the chiral perturbation results in agreement 
with that obtained from the KN scattering length. 




P/Po 



Fig. 3 Kaon and antikaon in-medium masses as functions of density from the mean-field approximation to the 
chiral Lagrangian. 

Antikaon properties in nuclear matter deserve some further discussions. The possibility of K~ condensation in 



neutron stars has been the focus of many recent studies 11C.146-14S]. The isospin averaged K N scattering length 



is negative in free space [ 15C ] , implying a repulsive K optical potential in the simple impulse approximation. However, 
a systematic analysis of the kaonic atom data shows that the K~ optical potential is deeply attractive, with a value of 



about 200± 20 MeV at normal nuclear matter density [ 151 1 . Unlike the kaon which interacts with nucleons relatively 
weakly so that the impulse approximation is reasonable, the antikaon interacts strongly with nucleons so we do not 
expect the impulse approximation to be reliable. The antikaon-nucleon (KN) interaction at low-energy is strongly 
affected by the A(1405) which is a quasi bound state of an antikaon and a nucleon in the isospin 1 = channel and 



can decay into T,tt channel [152 153 1. Thus, in principle one needs to carry out a coupled-channel calculation for 
KN, An and T,tt including the effects of A(1405) in both free space and in nuclear medium [152 |l57 |. Because of 
Pauli blocking effects on the intermediate states, the possibility of forming a bound KN state (A(1405)) decreases 
with increasing density, leading to a dissociation of A(1405) in nuclear medium. This induces a transition of the K~ 
potential from repulsion at very low densities to attraction at higher densities [ 153| , 156[ . It was shown in Ref. [157|, 
that, because of the diminishing role of A(1405) in dense matter, the prediction for the kaon condensation threshold 
is relatively robust with respect to different treatments of A(1405). Also, the strong attraction found in kaonic atoms 
can be understood in chiral perturbation theory if the possible scaling of the pion decay constant / in nuclear medium 
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is included @[u||l45) 



C. vector mesons 

Based on the restoration of scale invariance of QCD, Brown and Rho have argued that masses of nonstrange vector 



mesons would be reduced in dense matter 158 15£ ] . This is supported by studies based on the QCD sum- rule approach 



L3|, the effective hadronic model including NN vacuum polarization 160 162] ], and the quark-meson coupling model 



163 1 



In the QCD sum-rule study of vector meson masses in nuclear medium, one considers the correlation function of a 
vector current, 

IW?) = i J e i<lx (TJ^x)M0)) p d 4 x, (33) 
In nuclear medium, it can be written as 

n 00 = q 2 n L , 

rr fx W\n , QWj'^tt fiA\ 
U v = [Sv - -^r) n T + (34) 

where and II ^ denote the transverse and longitudinal parts of the polarization tensor, respectively. At zero 
momentum, q = 0, there is no distinction between the transverse and longitudinal directions, and Hl is thus related 
to ny, i.e., 

n T = (lu 2 - q 2 )n L , 

^L = -^9^U^ . (35) 

In this case, only the longitudinal correlation function is needed. 
For a rho meson, the current is taken to be 

jjp) = *r»«-*r» d . (36) 

For large Euclidean four momenta, Q 2 (= —q 2 = — s) — > oo, Hl(Q 2 ) can be evaluated perturbatively by the operator 
product expansion (OPE) [ p8| . Including operators up to dimension 6, one has 

n iW *) - ^ (i + .o g | + 'p m)e + 5^(0^0-), 

-^(({7,7 5 A"«)(57'75A''4)) P - ^ Y. Vft^'titi-f >•'</)),■ (37) 

q'—u,d 

In the above, a s is the QCD coupling constant while Qo is an arbitrary scale parameter. Both light quark masses 
and their expectation values are taken to be the same, i.e., m q = m u w rrid and (qq) p — (uu) p ~ (dd) p . 

The imaginary part of the correlation function at s > is related to the rho meson spectral function in medium 
and is normally parameterized phcnomcnologically by a contribution from the rho meson pole and a continuum, i.e., 

87rImn L (s) = FS(s - mf) + (l + — ) 9{s - «„), (38) 
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where F is a constant. 

The theoretical side (Q 2 > 0) and the phenomenological side (s > 0) are related through the dispersion relation, 

Ren L (Q 2 ) = - f ds ^— ^ + subtractions. (39) 

7T Jo S + Q 

From the above equation, the rho meson mass m* and the continuum threshold sq can then be determined. 

To suppress contributions from the higher order operators, one usually introduces the Borel transform defined by 

1 / rl \ n 

Lm b = lim (Q 2 )« - , (40) 

Q 2 ,rwoo;Q2/n=Mf, (t7-1)P \ dQ 2 ) 



where Mb is the Borel mass in Eq. ( |ll[) . The Borel transform also removes the need for subtractions in the dispersion 
relation. 

Carrying out the Borel transform of both sides of eq. ( |39| ) and taking the ratio of the resulting equation to its 
derivative with respect to — 1/Mg, one obtains 

m * (i + ^)[i - (i + ity- sa,Ml \ - *-itm P - n&riGrcnp + iSf « 



(41) 



In deriving eq. (^) , one has made use of the mean-field approximation, 



{{qi^ a q){qri^ a q)) P « -(qq)% 



16 

y 
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{{qi^ a q){q'ry i q'))p « (42) 

Also, terms involving nonscalar operators due to lack of Lorentz invariance in a medium are not shown. Similar sum 
rules can be derived for the omega and phi meson masses, using the currents — (u^ + d^y ^d) / '2 and Jp — S7 M s, 
respectively. 

To determine the rho meson mass, one minimizes the sum rule, Eq. (fll|), with respect to the Borel mass Mb at 
an optimal sq. In free space, using the commonly adopted values for the vacuum condensates ((qq)o = — (245 MeV) 3 , 
^(G /Jly G M ' y )o = (330 MeV) 4 ), quark mass (m q = 5.5 MeV), and QCD coupling constant (a s = 0.3), the empirical rho 
meson mass (776 MeV) can be reproduced with so ~ 1.77 GeV 2 in the range of Borel masses 0.55 < M B < 0.75 GeV 2 . 

In nuclear medium, vector meson masses depend on in-medium condensates. Because of small light quark masses, 
terms in Eq. (^l[) linear in quark condensates can be neglected. The results for p and ui mesons in nuclear medium thus 
depend on the density dependence of four-quark condensates. For phi meson, (ss) p dominates due to the relatively 
large strange quark mass m s ~ 100 ~ 200 MeV, its in-medium mass thus depends on the nucleon strangeness content 
as (ss)p ~ (ss)o + (N\ss\N)pn, where (ss) w —0.8(qq)o is the strange quark condensate in vacuum. Based on 
the factorization assumption for the four-quark condensate, i.e., Eq. (p2|), and using a simple delta- function plus 
continuum ansatz for the vector meson spectral function as shown in Eq. (|38|), Hatsuda and Lee have obtained the 
following results for the in- medium vector meson masses Jl3| , 

m* 



777 



1 - 0.18(pAr/p ), (43) 



and 



mi 

J «l- 0.15y( PN /p ). (44) 

777^, 
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The density dependence of rho (omega) and phi meson masses is shown in Fig. 4, where the nucleon strangeness 
content is taken to be y — 0.17. It is seen that the rho (omega) meson mass decreases significantly with density due 
to a strong density dependence of the light quark condensate. Since the strange quark condensate does not change 
much in nuclear medium as a result of the small nucleon strangeness content, the phi meson mass thus shows a weaker 
density dependence. However, the phi meson mass would decrease substantially if the temperature is high so that 



many strange particles are present and reduce thus the strange quark condensate [164|. The temperature dependence 
of phi meson mass is shown in Fig. 5 and is seen to decrease significantly at high temperatures. We note that 



recent studies on phi meson mass at finite temperature using the hidden gauge theory 165. 166] show a much smaller 
reduction than that from the QCD sum rules. This may be due to the incomplete treatment of vacuum effects in the 
latter approach as it is carried out only at the one-loop level. 



s~ aoo - 
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Fig. 4 Rho and phi meson masses as functions of density from QCD sum rules, (from Ref. fil 



A similar QCD sum-rule calculation of in-medium rho-meson mass has been carried out by Asakawa and Ko [167] 
using in the RHS of Eq. (|39| ) a more realistic spectral function that is obtained by coupling the rho meson to pions via 
the vector dominance model and using the pion in-medium dispersion relation determined from the delta-hole model 



as discussed in Section II.B.l [ 168| , 169[ . A decrease of rho meson mass with increasing density as that of Ref. [jl3| has 



also been obtained. This is due to the stronger effect from reduced quark condensate in the presence of nucleons than 
that from the matter polarization. 



Vector meson masses have also been studied using effective hadronic models 1 16C - 162 ] . Taking into account the 
change of NN vacuum polarization in medium as a result of reduced nucleon mass, the rho and omega meson masses 
have been found to decrease with density and temperature as in QCD sum-rule approach. If one includes only the 
matter effect, i.e., the polarization of the Fermi sea, the rho meson mass is seen to increase at high density. The 
inclusion of the vacuum effect, i.e., the polarization of the Dirac sea, however, brings the rho meson mass down in 
medium. The modification of the vacuum polarization in effective hadronic models is thus related to the change of 
quark condensate in medium in the QCD sum-rule approach. Since the nucleon is a composite particle, the large 
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effect of NN vacuum polarization has, however, been questioned in Ref. [ 17C ] . 
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Fig. 5 Phi meson mass as a function of temperature from QCD sum rules, (from Ref. |164|) 



Experimentally, the quenching of the longitudinal response in quasielastic electron- nucleus scattering [171], and 



the enhancement of the -ftT + -nucleus scattering cross section [172,173| have been considered as indirect indication of 
the decrease of omega meson mass in nuclear medium. More direct observation of the vector meson mass in dense 
matter can be provided by the dilepton invariant mass spectrum from heavy-ion collisions. Since low-mass dileptons 
are mainly produced from pion-pion annihilation that proceeds through a rho meson, a change in the rho-meson mass 
in dense matter is expected to be reflected in the dilepton spectrum as a shift of the rho-meson peak to a lower mass 
|58[ . Similarly, the in-medium properties of phi meson can be studied from the dilepton spectrum via kaon-antikaon 
annihilation. 



III. THE RELATIVISTIC TRANSPORT MODEL 



A. a brief review of the development of transport models 



Heavy-ion collisions involve very complicated nonequilibrium dynamics and exhibit different features at different 
incident energies. At low energies, because of Pauli blocking of two-body collisions, a suitable approach is the mean- 



field theory, such as the time-dependent Hartree-Fock (TDHF) theory and its semiclassical approximation [174-178]. 
At high energies, the reaction dynamics at the initial stage is mainly governed by two-body collisions, so the mean-field 
can be safely neglected. Cascade-type models, in which heavy-ion reactions are visualized as a sequence of independent 
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two-body collisions, have been developed for this purpose, first for heavy-ion collisions at Bevalac energies [179,180], 



and recently for heavy-ion collisions at AGS |181] and SPS energies 182 



At energies between these two limits, and for the expansion stage of heavy-ion collisions at high energies, both the 
mean field and two-body collisions are important and need to be included in the model. There have been attempts 
to extend the TDHF approach to include the effects of two-body collisions. In this way, one obtains the extended 
time-dependent Hartree-Fock (ETDHF) equation, which were, however, found to be too complicated to be useful 



for heavy-ion collisions 183,184]. On the other hand, semiclassical approximations, which make it possible to treat 



heavy- ion collisions more efficiently, have been shown to give results that are close to that of TDHF [185]. This has 
thus led to the development of the well-known Boltzmann-Uehling-Uhlcnbcck (BUU) equation, which is also known 
as the Vlasov-Uehling-Uhlenbeck (VUU), Landau- Vlasov (LV), and Boltzmann-Nordheim-Landau (BNL) equations 



[ 186 - 191 1 . The BUU equation includes simultaneously the effects of mean field, two-body collisions, and the Pauli 
principle. Its solution then gives the time evolution of the one-body phase space distribution function. 



The BUU equation is usually solved by the test-particle method 1 192 ] , in which the propagation of particles in the 
mean field (the Vlasov part) is given by the Hamilton's equation of motion, while two-body collisions (the Uehling- 
Uhlenbeck part) are treated by Monte-Carlo procedure. The introduction of test particles is necessary to remove 
numerical fluctuations associated with the evaluation of the mean field, but at the same time they also suppress 
correlations that are important for describing the instability and fragmentation of the system at low densities. A 
similar approach, the quantum molecular dynamics (QMD) (24 , 193 ,194] , has thus been developed. In the QMD model 
each particle is represented by a Gaussian wave packet in both space and momentum. In this way, one can obtain 
a smooth mean field without introducing test particles (in other words, one test particle for one physical particle). 
The QMD model seems to provide a reasonable description of the nuclear fragmentation phenomena p4[ . Similar 
but improved models, such as the fermionic molecular dynamics (FMD) ]195| and the antisymmetrized molecular 



dynamics (AMD) [196,197], have also been developed. 

The two important ingredients in these transport models, namely, the mean-field potential in the Hamilton's 
equation of motion and the in-medium cross sections for two-body scattering, are usually introduced separately. For 
example, the Skyrme parameterization |2llp4| for the mean field and the Cugnon parameterization for the NN cross 



sections [21 198 are often used in transport models. Also, the momentum dependence of the nuclear optical potential 
has been introduced as it has been found to be important in describing the nucleon collective flow in heavy-ion 



collisions [199 202]. Since the mean field and in- medium two-body cross sections are related through the in-medium 
effective NN interactions (the G matrix), using separate parameterization inevitably introduces ambiguities in the 
transport models. There have been attempts to derive consistently the mean field and the in- medium two-body cross 



sections from the same underlying NN interactions [ 203 - 205 1 . 

Also, relativistic effects are expected to become important in heavy-ion collisions at high energies. Although 
relativistic kinematics have been incorporated in most transport models, the relativistic covariance, or the frame 
independence, of the transport model has only been addressed in the relativistic quantum molecular dynamics (RQMD) 



[182 1 and similar models |206| by the use of Hamiltonian dynamics that is constrained by Poincare invariants. The 
meson-exchange nature of NN interactions, by which the nucleon mean-field potential can be separated into different 
Lorentz components (e.g., scalar and vector potentials), has been included in the relativistic Boltzmann/ Vlasov- 
Uehling-Uhlenbeck (RBUU/RVUU) approach ]2g|-|3"o[] . The relativistic transport model allows one to investigate 
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consistently the medium effects on hadron properties through the change of the scalar and vector potentials. We shall 
discuss briefly in the next subsection the derivation of the relativistic transport model from Walecka-type models. 



B. the relativistic transport model 

Based on the Lagrangian of nonlinear a-uj model, Eq. (^), in which the nuclear matter is treated as a system of 
interacting nucleons and mesons, a transport equation for the phase space distribution function of nucleons has been 
derived in Refs. p2]p5| , p8 , 207 - 210 ] using the mean-field approximation.. 



An important quantity in these derivations is the nucleon Green's function, defined by 

iG(x u xy) = (\T[N(xi)N(xy)]\) 

= 9(h -t v )iG > (x 1 ,x v ) + 0(t v -t 1 )iG < (x 1 ,x v ), (45) 

where (• • •) denotes the expectation value in the nuclear many-body state and T is the time-ordering operator defined 



on a contour in the complex time plane [211,212]. It satisfies the following equation 



(*7m^x! ~ m N )G(xi, x v ) = 5(xi - xy) + j d uE(xi, u)G(u, Xy), (46) 

where £ is the nucleon self-energy and can be evaluated via the perturbative expansion. The first-order nucleon 
self-energy is normally written as Yi{x\,X2) ~ ^hf(xi)8(xi — X2) with 



{x) = V s +i"{^ v {x). (47) 



The scalar self-energy can be absorbed into the nucleon effective mass as in Eq. (|J) while the vector self-energy leads 
to an effective momentum p* = + Xv>. In this limit, the nucleon Green's function satisfies the equation, 

- " m* N (x 1 )}G < (x 1 ,xy) = 0. (48) 

Introducing the Fourier transform of the nucleon Green's function 

G<(x,p*)= f d i ye^G<(x + y/2,x~y/2), (49) 



then in local density and semi-classical approximations, it is related to the seven dimensional on-shcll nucleon's 
phase-space distribution function /(x,p*), i.e., 

Tr[G<(x lP *)] = ^167r5(p* 2 - m* N 2 )m* N f{x,p*). (50) 

From Eq. p^), a relativistic Vlasov equation follows and has the form 

{[d£ - (d»K - d^ v )df] P ; + mUd^m* N )df}f(x, p*) = 0, (51) 



This equation can be more conveniently written as [213] 

d 



dt 

where 



/ + v-V x /-V s C/-V f ,/ = 0, (52) 
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v = p*/E*, U = E* + {g^/m^fpN, (53) 

with E* = (p* 2 + m*' 2 ) 1 / 2 . 

In terms of the phase-space distribution function, the scalar and current densities in the local-density approximation 
can be expressed, respectively, as 

Ps=4 J §^K x >P*> l */ E *> P ^ 4 J S^K X >V*K/ E *- ( 54 ) 

In heavy-ion collisions the nucleon vector potential has both spatial-like and time-like components and are related to 
the nuclear current density p M . 



The relativistic transport equation is solved by the test-particle method [ 192 1 , so each nucleon is replaced by a 
collection of test particles. The one-body phase-space distribution function is then given by the distribution of these 
test particles in phase space. To solve the Vlasov equation is thus equivalent to solving the following classical equations 
of motion for all test particles, 

f| = p*/E\ ^ = -V x U(x). (55) 



These equations of motion can be extended to other hadrons. For baryon resonances, they are the same as Eq. (|5J 
as we have assumed that they have same mean-field potentials. For lambda and sigma hyperons, the corresponding 
mean- field potentials are taken to be 2 /3 of those for the nucleon. The equations of motion for an antiproton are the 
same as those of the nucleon except that the sign of the vector potential is changed due to G parity. For a kaon, 
the equations of motion in the mean-field approximation to chiral Lagrangian p6[ are obtained from Eq. (pq) with 
E* = m 2 K + p 2 ^jr-ps + (| jt) and U = to(p, pb) — ^o(p), where uj(p, pb) and w (p) are the kaon dispersion 



relation in medium and in free space, respectively [ 143 1 . Again, the vector potential for an antikaon has an opposite 



sign from that of the kaon 49 



If one includes the second-order Born terms in the nucleon self-energy, then its imaginary part can be shown to 
lead to a collisional integral on the right hand side of the relativistic Vlasov equation p8| , |208|] , i.e., 

/*■&■<*• +ps-ps-po 

•{/(*, P3)/(*> pS)[i - /(*> p*)][l - f(x, p5)] - f(x, p*)f(x, p5)[l - f{x, p%)] 

■[l-fix,pt)}}, (56) 

where v is the relative velocity between the colliding particles and 4£ is the two-body NN differential scattering cross 
section calculated from meson exchanges. The factor (1 — /) takes into account the Pauli-blocking of fermions in the 
final state of a scattering. 

For heavy-ion collisions at a few GeV/nucleon, nucleons, delta resonances, and pions are the most important degrees 



of freedom. As in normal BUU model the isospin-averaged cross sections measured in free space [ 198 are used 



for the elastic (NN — ► NN) and the delta excitation (NN — > NA) process. The cross section for the inverse process 



7VA — > NN is determined from the detailed balance relation [ 214 1. Both the nucleon-delta (7VA — > NA) and the 
delta-delta (AA — > AA ) elastic collision are also allowed, and their cross sections are assumed to be the same as 
that for the nucleon-nucleon elastic scattering at the same center-of-mass energy. 

When a delta is formed, its mass distribution is taken to be of the Breit-Wigner form 
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P(mA) -(™ A -mo)2 + (r (g )/2) 2 ' (5?) 



where too = 1-232 GeV and the momentum-dependent delta width [215] , r(g), is 

0.47g 3 

[i + o.%7^) 2 ]to2 • 



r (g)= M ■ n^l E ~ (58) 



In the above, q is the pion momentum in the rest frame of a delta and is related to its mass toa by q — 
J [to a - (m N + to^) 2 ] [m\ - (m N - m^) 2 ] /(2to a ). 

The collision between two baryons is treated in the same way as in cascade model [198|. A collision occurs when the 
distance between them is less than \J a /n with a being their interaction cross section. After the collision, directions of 
their momenta change in a statistical way according to the angular distribution. However, collisions are allowed only 
among particles in the same simulation but the mean nuclear density and current are computed with test particles 
from all simulations in the ensemble. 



We note that a different treatment of two-body collisions has recently been introduced in Ref. 1 216 ] . In this approach, 
the collision between two particles is treated classically so that not only the angular momentum is conserved but also 
the reaction plane is preserved. Furthermore, only repulsive orbits arc allowed in the collision due to the largely 
repulsive force in nucleon-nucleon scattering below 1 GeV/nucleon. An attempt to address the question of separating 
the nuclear interaction into a mean-field and a two-body collision part in nuclear collisions has recently been given in 



Ref. [217 1 . 

Pion production is included in the relativistic transport model through delta decay, i.e., A — > ttN, using the 
momentum-dependent delta decay width. The inverse reaction ttN — ► A is also included to take into account pion 
absorption [218- 22C[ . In most transport models for pions, one does not include the pion potential so they propagate 
as free particles in nuclear medium. However, medium effects on pions are included in Ref. |53| by treating them 
as quasiparticles. It is found that softening of the pion dispersion relation might be responsible for the observed 



enhancement of low transverse energy pions in heavy-ion collisions at both Bevalac and SIS energies |)4 55 . 

Because of their small production probability in hadron-hadron collisions, other particles, such as the kaon, antikaon, 
antiproton, and dilepton, are treated perturbatively, i.e., the collision dynamics is not affected by their presence 
p2l] , p22| . These particles may also scatter with other hadrons. A kaon interacts only elastically with a baryon 
due to strangeness conservation, and the kaon-nucleon scattering cross section at low energies is about 10 mb in 



free space [223]. On the other hand, an antikaon, in addition to clastic scattering, can be absorbed by a nucleon 



through the strangeness-exchange process KN — ► irY. These cross sections are also taken from Ref. [223]. Similarly, 
antibaryons can be annihilated by baryons, and their cross sections are taken from Ref. p24| . To include the final- 



state interactions of these particles, a perturbative test particle method has been introduced in |225|. In this method, 
one allows many particles to be produced in a hadron-hadron interaction. Each produced particle is then assigned 
a production probability, which is given by the ratio of its production cross section to the hadron-hadron total cross 
section. Their motions are then followed by solving the classical equations of motion. The collisions of these particles 
with nucleons are treated as usual except that their effects on nucleons are neglected, i.e., one does not allow the 
nucleon momenta to change in such collisions. 

With reduced hadron masses in medium, the above cross sections should be modified as well. For nucleon-nucleon 



elastic scattering, the in-medium cross section has been evaluated in Ref. [226 227] in the DBHF approach using the 
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one-boson-exchange potentials. It is found that the magnitude of the cross section decreases with density and the 



differential cross section becomes more isotropic in medium. Experiments on proton- nucleus scattering |228 230] seem 
to be consistent with these predictions. The NN — ► NA cross section in medium has been studied in Refs. 231 - 233 
via the one-pion exchange model. Including the softening of the pion dispersion relation, its magnitude is found to 
increase with density. Recently, the same model has been used for evaluating the AA — > AA cross section in medium 
p34j . Although the magnitude of the cross section is found to decrease with density, the density dependence is less 
pronounced than that of the nucleon-nucleon elastic cross section. 

The A «-> ttN cross section in medium can be easily evaluated using in-medium hadron masses if one assumes 
that the cross section remains to be of the Breit-Wigner form. No theoretical studies have been carried out for the 
kaon-nucleon clastic cross section in medium. It is thus of interest to pursue such a study following the method of 



Rcf. [227] by using the meson-exchange model of Ref. 1 235 ] . Similar calculations should also be done for the antikaon. 



The production cross sections for kaon, antikaon, antiproton, and dilepton from hadron-hadron interactions in both 
free space and medium will be discussed below. 

IV. MEDIUM EFFECTS IN HEAVY-ION COLLISIONS 

In this section, results obtained from the relativistic transport model are presented and compared with the exper- 
imental data from heavy-ion collisions. In particular, discussions will be given for particle production and flow at 
SIS/GSI, strangeness enhancement at AGS/BNL and SPS/CERN, and dilepton production in heavy- ion collisions. 

A. heavy-ion collisions at SIS energies 

In heavy-ion collisions at SIS energies (~ 1 GeV/nucleon), the colliding system consists mainly of nucleons, delta 
resonances, and pions. In addition, particles with small production probability, such as the eta, kaon, antikaon, 
antiproton, and dilepton, can be created from baryon-baryon, meson-baryon, and meson-meson interactions during 
the course of the collisions. Since the threshold energies for eta, kaon, antikaon, and antiproton production in NN 
collisions in free space are, respectively, 1.26, 1.58, 2.5, and 5.6 GeV, their production in heavy-ion collisions at SIS 
energies is mostly below the production threshold in NN interactions in free space, and is known as 'subthreshold 
particle production'. Subthreshold particle production offers the possibility of exploring the extreme region of phase 
space, and thus may carry valuable information about the early, violent stage of heavy-ion collisions during which 
both density and temperature are high. In particular, the mean-field potential has been found to play an important 
role in subthreshold particle production p5jppg|j5l],|5l, p36l - p4^| . This can be seen from Fig. 6, where we show the 



density dependence of the Q-value, which is defined by the difference between total hadron masses and potentials in 
the final and initial states, for the reaction NN -> NAK, NN -> NNKK, and NN -> NNNN. They are evaluated 
with the in-medium hadron masses and potentials discussed in Section II for the soft nuclear equation of state in 
Table I. It is seen that the Q-value increases slightly for the reaction NN — > NAK but decreases appreciably for the 
reactions NN — ► NNKK and AW — > NNNN. In the ideal case that the scalar and vector potentials on a hadron 
in nuclear medium is determined by the sum of scalar and vector potentials on its constituent quarks, the Q-value for 
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the reaction NN — > NAK is independent of density, while it is 2/3£s and 2Sg, respectively, below the values in free 
space for the reactions NN — > NNKK and NN — > NNNN . At twice normal nuclear density, their values are -177 
MeV and -530 MeV, respectively, which arc slightly less than that shown in Fig. 6. 
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Fig. 6 The density dependence of the Q-value for the reactions NN — > NAK (solid curve), NN — > NNKK 
(dashed curve), and AW _> NNNN (dotted curve). 



i. subthreshold kaon production 



The first experiment on kaon production in heavy-ion collisions was carried out at Bevalac in the eighties at incident 



energies around 2 GeV/nucleon [ )244 ,245|. These experiments had generated many theoretical studies based on both 
cascade and transport models |200p2l|p2^p3^p4^ - p4^ ] . Recently, there are new experiments on kaon production in 
heavy- ion collisions at incident energies around 1 GeV/nucleon by the KaoS collaboration at GSI p8|| . This has led 
to a resurgence of theoretical studies based on both relativistic and nonrelativistic transport models 



17|,|237H240|. 

The study of kaon production in heavy-ion collisions at subthreshold energies allows us to investigate the properties 
of the dense matter formed in the initial stage of the collisions. These include the nuclear equation of state, the role 
of baryon resonances, and the kaon in- medium properties. 

To study the kaon production cross section in heavy-ion collisions, we need its production cross section in baryon- 
baryon interactions as the contribution from other processes is small. For example, it has been shown that in heavy-ion 



collisions the contribution to kaon production from pion-baryon interactions is only about 25% 24E ] while that due 



to many-body collisions is about 10% |250| . Randrup and Ko [221] have analyzed the available experimental data 
and proposed the following parameterization for the isospin-averaged kaon production cross section from the NN 
interaction in free space 



C r NN^BYK+ 



{y/s) M 36 



Pn 



m K 



fib, 



(59) 
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where the kaon maximum momentum p max is related to the NN center-of-mass energy y/s by p max = 
yj [s — (m B + m Y + m K ) 2 ] [s — (m B + m Y — m K ) 2 ]/As. 

Kaon production cross sections from the nucleon-delta and delta-delta interactions have also been analyzed in Ref. 
J221[ based mainly on isospin arguments. It has been found that the following scaling relations approximately hold, 

3 

< J NA^BYK+ (vS) ~ -^Cr NN ^ BY K+ (Vf)j 

&AA^BYK+(Vs) ~ ^NN^BYK+iVs)- (60) 

In addition to the total kaon production cross section, one also needs the kaon momentum distribution from baryon- 



baryon interaction in free space. This has been parameterized in Ref. 221 | based on phase space considerations, i.e., 

da r ,12 / p \ 2 ( p 



dp Pmax VPmax/ \ Pn 



which reproduces reasonably the experimental data J221| . 

Using the one-pion plus one-kaon exchange model [251|, it has been recently shown that near kaon production 
threshold, where no experimental data are available, the cross sections in the Randrup-Ko parameterization and scaling 
ansatz are overestimated for the nucleon-nucleon interaction but underestimated for the nucleon-delta interaction. 
However, these discrepancy may be due to the neglect of final-state interactions, which have been shown to be 



important near the kaon production threshold [251]. 




Fig. 7 Kaon momentum spectra from Au+Au collisions at 1 GeV/nucleon, obtained with the soft (solid) and stiff 
EOS (dashed), respectively. The left, middle, and right panels are from Refs. Ji^| , [239], and [ 238 1 , respectively. The 
experimental data from the KaoS collaboration El] are also shown. 
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As first pointed out by Aichelin and Ko [ 222 ] , the kaon yield from heavy-ion collisions at incident energies below the 
kaon production threshold in NN interaction in free space is sensitive to the nuclear equation of state at high densities. 
For Au+Au collisions at 1 GeV/nucleon, this dependence is shown in Fig. 7 in which the kaon momentum spectra 
obtained in three different calculations using the soft and the stiff equation of state are compared with experimental 
data from the KaoS collaboration . The results in the left panel of Fig. 7 arc from Rcf . 46 based on the rclativistic 
transport model including medium modifications of the kaon properties. In this calculation, the kaon production cross 
section in medium is obtained by using in- medium hadron masses to evaluate the p max in Eq. (59). The results in 



the middle panel are from Ref. [ 239 1 based also on a relativistic transport model. The results in the right panel are 



from Ref. [ 238 ] using the nonrelativistic QMD model. Theoretical results from relativistic transport models with the 
soft EOS are in reasonable agreement with the experimental data, while that with the stiff EOS are below the data 
by about a factor of two. However, the results from the QMD calculation seem to favor a stiff nuclear equation of 



state. As explained below, the QMD calculation of Ref. [23S] has overestimated the contribution from delta particle 
as they are only allowed to decay at freeze out. Thus, the KaoS data from Au+Au collisions at 1 GeV/nucleon seem 
to be consistent with a soft EOS. 




Fig. 8 Kaon momentum spectra from Au+Au collisions at 1 GeV/nucleon, obtained with (solid) and without 
(dashed) the kaon scalar potential, (from Ref. 

Subthreshold kaon production in heavy-ion collisions is also sensitive to the kaon scalar potential, and this is 
illustrated in Fig. 8. The solid and dashed curves correspond to results with and without kaon scalar potential, 
respectively. It is seen that without the attractive kaon scalar potential, the theoretical results are about a factor of 
3-4 below the experimental data. This can be easily understood from Fig. 6 in terms of the Q-value for the reaction 
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NN — > NAK in medium. At nuclear density p = 2po, where most kaons are produced in heavy- ion collisions, the 
Q-value is about 33 MeV above that in free space. Without the attractive kaon scalar potential, which is about 90 
MeV, the Q-value is increased by the same amount. This increases the kaon production threshold and thus reduces 
its yield. 



In the calculation of Ref. [23£] based on the relativistic transport model, medium effects on kaons are not included 
so kaons are treated as free particles as in the nonrelativistic QMD calculation [ 237 , 238|] . The results also agree with 
the experimental data and are thus comparable to that of Ref. W5[, as shown in the middle panel of Fig. 7. Although 



the kaon mean-field potential is taken to be zero in Ref. [239], the scalar and vector potentials for the hyperon are 
taken to have the same strength as that for a nucleon, which is stronger than the hyperon mean-field potential used 
in Ref. Ea|, i.e., 2/3 of the nucleon mean-field potential. Thus, as far as the kaon production threshold is concerned, 



the treatment of Ref. [239] is approximately equivalent to Ref. []45[ in which both kaon scalar and vector potentials 
are included. 




Fig. 9 Contributions to the total kaon production probability from the NN, NA and AA interactions in Au+Au 
collisions at 1 GeV/nucleon and impact parameter b=0 fm. (from Ref. p5[) 

Most kaons in heavy-ion collisions at subthreshold energies are produced from the two-step process in which a 
nucleon is first excited to a resonance and a kaon is then produced from the interaction of the resonance with another 
baryon. At incident energies around 1 GeV/nucleon, only the delta resonance is appreciably produced. Subthreshold 
kaon production can thus serve as a good probe of the properties of (delta) resonance matter. In Fig. 9, the total 
kaon production probability in Au+Au collisions at an incident energy of 1 GeV/nucleon and an impact parameter of 
b=0 fm from the relativistic transport model is separated into contributions from the nucleon-nucleon, the nucleon- 



2G 



delta and the delta-delta interaction. It is seen that the contribution from the nucleon-delta interaction accounts 
for more than half of the total kaon yield, and the sum of contributions from the nucleon-delta and delta-delta 
interactions is about 75% of the total kaon yield. Contrary to that of Refs. [|237 , |238 |, the contribution from the 
delta-delta interaction is slightly smaller than that from the nucleon-nucleon interaction. This is mainly due to 



following two reasons. First, in Refs. [237, 238|, deltas are allowed to decay only at final stage of the collision. This 
treatment certainly overestimates the contribution from the nucleon-delta and delta-delta interactions. Secondly, 



kaon production in Refs. [237,238] is calculated in the QMD with a Skyrme-type momentum-independent mean-field 
potential. It is well-known that introducing a momentum-dependence in the nucleon mean-field potential, which is 



automatically included in the relativistic transport model [252|, reduces the number of two-body collisions and hence 
the delta abundance, thus leading to a smaller contribution from the delta-delta interaction as compared to the case 
with a momentum-independent mean-field potential. 

Although the interaction of a kaon with a nucleon is relatively weak as compared to other hadrons, they are still 



scattered by baryons in the dense medium [225,253]. Furthermore, kaons propagate in both scalar and vector mean- 
field potentials. These final-state interactions change the kaon momentum spectra but not its total yield, leading to 
an increase of both its yield and the slope of its spectra at large angles. Actually, final-state interactions help to bring 
the theoretical results at #i a b = 44° in better agreement with the experimental data. Since final-state interactions 
increase the kaon yield at large angles, it would be of interest to have experimental data to test this prediction. 



2. subthreshold antikaon production 

Effects of the attractive scalar potential can be more clearly seen in antikaon production than in kaon production 
as it leads to a significant reduction of the Q- value for the reaction NN — > N N K + K~ in nuclear medium as shown in 



Fig. 6. From the analysis of available experimental data, Zwermann and Schurmann |254| have proposed the following 



parameterization for the antikaon production cross section from nucleon-nucleon interaction in free space 

dp VPmax/ V P 



d ° ~ 0.75 (-*-)' ( 1 - -E-Y m b/GeV, (62) 



where p max = y[s — (2m + 2m#-) 2 ](.s — Am 2 ) / (As) . It has been usually assumed that the same cross section is 
applicable for antikaon production from the delta-nucleon and the delta-delta interaction. 

Including all final-state interactions (propagation in mean fields, elastic scattering, and absorption) for antikaons, 
the results from the relativistic transport model for the antikaon production cross section from 58 Ni+ 59 Ni collisions at 
1.85 GeV/nucleon are shown in Fig. 10 j|9|. The dashed curve gives the results using free kaon and antikaon masses, 
while the results with in-medium kaon and antikaon masses are shown by the solid curve. Also shown in the figure 
by solid circles are the experimental data from the SIS at GSI |5(| . 

Using free kaon and antikaon masses, the theoretical results are about a factor of 5-10 smaller than the experimental 
values. The antikaon yield is enhanced by a factor of 4-5 when in-medium kaon and antikaon masses are used. The 
theoretical results are now in reasonable agreement with the experimental data, except at p^ 6 =1.0 GeV/c where they 
are still below the data by about a factor of two. Including possible decrease of the pion decay constant in medium, 
which leads to a stronger reduction in antikaon mass, is expected to further increase the calculated antikaon yield. 
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Subthreshold antikaon production in heavy-ion collisions has also been measured earlier at Bevalac for Si+Si 
collisions at 1.55 and 2.1 GeV/nucleon |245|. Theoretical calculations using the QMD model without medium effects 



give results that are also smaller than the experimental data by about a factor of 3-4 [24£ ,255]. Again, this discrepancy 
is likely explained by the modifications of kaon and antikaon masses in medium. 




Fig. 10 Antikaon momentum spectra from Ni+Ni collisions at 1.85 GeV/nucleon, obtained with free (dashed) and 
in-medium (solid) kaon and antikaon masses. The experimental data from Ref. |5(j|] are shown by solid circles, (from 
Ref. @) 

It should be noted that an antikaon can also be produced from the process ttY — > KN between a pion and a 



hyperon Y produced from the reaction BB — > NYK. However, with the cross section of Ref. 1 256 obtained from the 



K-matrix analysis of Ref. [223], it has been found that the antikaon yield from ttY — > KN in heavy-ion collisions is 



only about 25% of that from BB — ► NNKK. This is less than that from calculations based on the kinetic model [257] 
and the cascade model |258f , where more than 50% of the antikaons are from this process. These earlier calculations 
are, however, less accurate than the one based on the transport model as they do not treat properly the collision 
dynamics and antikaon absorption. Similarly, antikaon production from -kB — > NKK has been found to be much 
smaller than that from BB interactions. Furthermore, the reaction NY — > NNKK is unimportant due to the small 
cross section [ |256t . We note that these estimates were all based on calculations without including medium effects and 
need to be further investigated. 

Antikaon production from decay could also be included in the transport model calculation. Unfortunately, the 
elementary 4> production cross section in nucleon-nucleon interaction is only scarcely available at high beam momenta. 
At a beam momentum of 10.0 GeV/c, the 4> production cross section is 1.0±0.1 jib [E59J in proton-proton interaction. 
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This is to be compared with the antikaon production cross section of 33±16 fib at the same beam momentum [ 260 1 . 
Also, recent experiments at the AGS indicate that, at an incident energy of 14.6 GeV/c per nucleon, the ratio of <j> to 
K ~ in a Si+Au collision is about 12% [(26l[ . We expect an even smaller ratio in heavy-ion collisions at beam energies 
of about 1-2 GeV/nucleon, and the contribution from phi decay to subthreshold antikaon production can thus be 



neglected 262 1. 



3. subthreshold antiproton production 



With a threshold at 5.6 GeV in nucleon- nucleon interaction in free space, antiproton production in nucleus- nucleus 
collisions at energies of a few GeV/nucleon is clearly the most extreme subthreshold process in particle production 
from heavy-ion collisions. The first observation of subthreshold antiproton production in proton-nucleus collisions 



dated back to the fifties ]263| ] and sixties [ 264, 265 1. The experiments at the Bevalac 1 245 and the JINR [266] in the 
eighties provided the first evidence of subthreshold antiproton production in nucleus- nucleus collisions. Recently, new 



measurements of subthreshold antiproton production have been carried out at both KEK [ 267 ] and GSI |50| . 

The antiproton momentum spectrum in proton-proton interaction in free space has been parameterized as 
p63| , p64| , p68 | 

2 / \ 2 
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(63) 



where the maximum momentum of the antiproton, p max , is related to the available center-of-mass energy y/s of the 
proton pair by p max = y/(s- 16m 2 N )(s - Am 2 N )/As. 

This parameterization is based on phase space arguments and describes reasonably well the measured antiproton 



momentum spectrum in proton- proton collisions [263,264 268 1. The total antiproton production cross section in Eq. 
(p3|) is fitted to the available experimental data, i.e., 



(TptVs) ~ 0.012(Vs - 4m A r) 1 - 846 mb. 



(64) 



The first investigation of subthreshold antiproton production using the transport model has been carried out by 



Batko et al. [26E]. In this study, the A(1232) degree of freedom is included and is found to play an important role. 
Neglecting the antiproton mean field and assuming that 90% of the primordial antiprotons are absorbed, Batko et al. 
[ p69| have achieved a reasonable description of antiproton production in proton-nucleus collisions, but their results 
for nucleus-nucleus collisions are about a factor of five too small compared with the experimental data. A similar 



calculation has been carried out by Huang et al. [ 255 1 using the QMD. Again, the antiproton yield in nucleus-nucleus 



collisions is less than the experimental data by about the same factor as in Ref. [26E]. In both calculations [ 255,269), 
antiproton annihilation is treated schematically, and the propagation as well as elastic scattering of antiprotons in 
medium are neglected. 
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Fig. 11 Antiproton momentum spectra from Ni+Ni collisions at 1.85 GeV/nucleon, and C+Cu collisions at 3.85 



GeV/nucleon. The left, middle, and right panels arc from Refs. p2|, [241], and 1243] , respectively. The experimental 



data from Ref. [pOj for Ni+Ni collisions and from Ref. [266 for C+Cu collisions are shown by solid circles. 



The importance of medium effects can be seen from the significant decrease of the Q-value shown in Fig. 6 for 
the reaction NN — > NNNN in nuclear medium as a result of the attractive scalar potential. This effect has been 
included in a number of studies based on transport models. In Fig. 11, theoretical results from these calculations 
for the antiproton differential cross section in Ni+Ni collisions at 1.85 GeV/nucleon and C+Cu collisions at 3.65 



GeV/nucleon are compared with the experimental data from both GSI [|50| and Dubna 266 . In the left panel, the 
results from Ref. |52j ] based on the relativistic transport model are shown. The antiproton mean-field potential is 
obtained from the nucleon one by G parity and is thus very attractive. It is indeed seen that the experimental data 
can only be accounted for if in-medium baryon masses (solid curve) are used. The results with dropping only the 
nucleon mass but not the antiproton mass are shown by the dashed curve and are a factor of 8 smaller. When free 
baryon masses are used in the calculation, the antiproton yield (dotted curve) is further reduced by a factor of 12. 
Overall, the antiproton production cross section is enhanced by about two orders of magnitude due to the reduction 
of baryon masses in medium. Similar results [5l| l have been obtained for antiproton production from Si+Si at 2.1 
GeV/nucleon ||| . 



The results in the middle panel of Fig. 11 are from Ref. [ 241 1 based also on the relativistic transport model. The 
results shown by the solid curve for the case with an in-medium antiproton mass are in good agreement with the 
data, while that with a bare antiproton mass (but still an in-medium nucleon mass) are below the data by about a 



factor of 5-10. We note that the theoretical results of Refs. [^2,241| agree with each other within a factor of two, and 



30 



the antiproton potential at normal nuclear matter density extracted from these calculations is in the range of -150 to 
-250 MeV. 



The results in the right panel of Fig. 11 are from Ref. 242 based on the nonrelativistic QMD model. Medium 
effects on the produced nucleon and antinucleon masses are taken from the Nambu— Jona-Lasinio model. It is again 
seen that without medium effects the theoretical results are much below the experimental data. These studies thus 
show that in order to describe the antiproton data from heavy-ion collisions at subthreshold energies it is necessary 
to include the reduction of both nucleon and antinucleon masses in nuclear medium. Even at AGS energies, which 
are above the antiproton production threshold in NN interaction, a recent calculation using the RQMD shows that 
medium modifications of the antiproton properties are important for a quantitative description of the experimental 



data (270 1 . 

On the other hand, there have been suggestions based on schematic considerations that the contributions from 



quasi-coherent multi-particle collisions [265] and meson-meson interactions (e.g. pp — > pp) [271] might be significant. 
More accurate estimates of these contributions using the transport model have not been done and are needed. 



4- soft pions 



Transport models have also been used to study pion production from heavy- ion collisions. Of particular interest 
is the experimental observation of enhanced low (transverse) momentum pions in these collisions at various energies 
[^4,55 272, 272]. Different mechanisms have been proposed to explain this phenomenon, including the collective flow 
|274| , the resonance decay at freeze-out [22C, 275- 277 , the finite pion chemical potential at freeze-out | 27ij , 279|] , and 



the medium modifications of the pion dispersion relation |53 28C-282]. In this subsection, we discuss mainly the 
results for heavy-ion collisions at Bevalac and SIS energies. 

Several nonrelativistic BUU calculations have been carried out by different groups for the pion spectra from heavy- 
ion collisions at 1-2 GeV/nucleon ]53] , 220] , 280 , 283 1 . The results of Ref. [|53| for the pion kinetic energy spectra from 
La+La collisions at 1.35 GeV/nucleon are shown in the left panel of Fig. 12. The dashed histogram is obtained 
using the free-space pion dispersion relation, and for low energy pions it differs significantly from the experimental 
data shown by solid circles with error bars j54j ]. With a softened pion in-medium dispersion relation described in 
Section II.B.l., there is an enhancement of low energy pions as shown by the solid histogram. This enhancement is 
due to several mechanisms. From delta decays in dense medium, quasipions in the pion branch generally carry less 
energy than free pions. With their increased widths, it is more probable to form small mass deltas, so that more 
soft pions are created from the decay of these deltas. Also, as quasipions propagate from a more dense region to 
a less dense region, as most likely happened during heavy-ion collisions, their momenta should decrease since the 
p-wave interaction reduces. Enhanced soft pions are also seen in more recent experiments at SIS |5q|. Calculations 



within transport models with a free pion dispersion relation fail to explain this enhancement [ 280, 283 1. A similar 



enhancement of low energy etas has also been observed [284|, and this may again be due to the change of the eta 
dispersion relation in medium as a result of the iV*(1535)-hole polarization. 
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Fig. 12 Pion transverse kinetic energy spectra in La+La collisions at 1.35 GeV/nucleon and in Ar+Ca collisions at 



1 GeV/nucleon. The left, middle, and right panels are from Refs. |53[, |220|, and |28C| ] , respectively. The experimental 
data from Ref. pj] for La+La collisions and from Ref. (55| for Ar+Ca collisions are shown by solid circles. 



It should be mentioned that there are still some controversies regarding transport model calculations of pion spectra 



in heavy- ion collisions. For example, in Ref. [ 220 enhanced soft pions can also be obtained by treating pions as free 
particles. This is shown in the middle panel of Fig. 12 for La+La collisions at 1.35 GeV. The low-momentum pions 
were found to come mainly from the decay of delta resonances at freeze-out. However, this was not confirmed in 



similar transport model calculations by other groups [53 280 283]. 

Furthermore, the effects of a modified pion dispersion relation in medium on pion spectra do not all agree. In Ref. 
|{280| , this effect is found to be very small due to a different treatment of the in- medium pion dispersion relation in 
the transport model. This is shown in the right panel of Fig. 12 where the solid and dashed curves are theoretical 
results obtained with the in-medium and the free-space pion dispersion relation, respectively. Therefore, further work 
on quasipions in transport models are required, and this has recently been taken up in Ref. ]285[ . 

Medium effects on pion (also kaon) momentum spectra in hot hadronic matter have been investigated in Refs. 
||281 , [282 | based on the 'optical potential' constructed from the forward scattering amplitude. The attractive pion 
collective potential was proposed in Ref. |28l|[ as a possible explanation for the observed enhancement of low transverse 



momentum pions (the soft pion puzzle) in heavy- ion collisions at SPS energies [273], in much the same way the 
attractive pion potential leads to enhanced soft pions in heavy-ion collisions at SIS energies. However, it was shown in 



Ref. [286 1 that in order for this attractive potential to be effective in cooling down pions, the colliding system has to 
expand slower than the speed of these low-momentum pions. Since at SPS energies, the system and hence the source 
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of the pion potential is made of pions as well, it is unlikely that the above condition can be satisfied. 



5. collective flows 



Besides particle spectra, the collective motion of particles and fragments, both in and out-of the reaction plane, has 
been extensively studied in heavy-ion collisions at SIS energies I|^,p7|, p0l| , p02p87| - p9^ | . It has been shown that the 
proton (and fragment) flow is sensitive to both the density (related to the nuclear equation of state) and momentum 



dependence of the nuclear potential [ 201 1 , as well as the in- medium NN cross sections [295- 297]. Pion flow has also 
been identified in heavy-ion collisions, and is seen to undergo a transition from flow at small impact parameters to 



antiflow at large impact parameters as a result of the shadowing effects from the spectator nucleons [283,292,293], 



Antiproton and antikaon flows in heavy-ion collisions have been seen in both the RQMD [ 294 1 and the relativistic 



hadronic cascade (ARC) [216] calculations. Because of their large annihilation cross section, both antiproton and 
antikaon Hows are found to be opposite to that of nucleons. Kaon flow has recently been proposed as another 
observable for studying the kaon properties in dense matter J5q], and is being studied experimentally by both the 



FOPI 1 57 1 and the EOS [28E] collaboration in heavy-ion collisions at a few GeV/nucleon. 




Fig. 13 Proton average in-plane transverse momentum from Ni+Ni collisions at 1.93 GeV/nucleon. The experi- 
mental data from the FOPI collaboration |^7j are shown by circles. 

In-plane particle flow is usually shown as the average in-plane transverse momentum (p x ) as a function of rapidity 
y. The results for proton flow in Ni+Ni collisions at 1.93 GeV/nucleon using the relativistic transport model with 
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a soft EOS are shown in Fig. 13. They are obtained after taking the impact-parameter-weighted average over b< 4 
fm, corresponding approximately to the centrality selection in experiments by the FOPI collaboration [^7| . Both the 
experimental data, shown in the figure by solid circles, and the theoretical results include a low transverse momentum 
cut of pt > O.brriN. It is seen that the theoretical predictions are in good agreement with the data except at the 
projectile and target rapidities where there are some differences, which are not understood at present. 

To extract nuclear equation of state at high densities, proton flow has also been studied extensively in nonrelativistic 
transport models with the inclusion of a momentum-dependence in the nuclear potential [201 202). The results of 



Refs. [201] and |202j arc shown in the left and right panels of Fig. 14, respectively. The flow parameter F in these 
work is defined by the average transverse velocity (p x )/m at y = 0. The momentum-dependence in the nuclear mean 
field was found to be important in reproducing correctly the dependence of the flow parameter on the incident energy 
and the colliding system. It is seen that the soft EOS with a compression modulus of about 200 MeV gives a good fit 
to the experimental data from both Nb+Nb and Au+Au collisions. 

In Fig. 15, pion flow in Au+Au collisions at 1 GeV/nucleon and impact parameters b =3, 6, 9 fm is shown. In 
central collisions (dotted curve), pions are seen to follow the flow of nucleons, while in mid-central (dashed curve) and 
peripheral (solid curve) collisions, the pion flow direction is opposite to that of nucleons. This transition from flow in 



central collisions to antiflow in mid-central and peripheral collisions has also been predicted in Refs. [283 ,293]. The 
'anticorrelation' of pions to nucleons at large impact parameters is mainly due to the rescattering and absorption of 
pions by the spectator nucleons. 




CO o 



Fig. 14 Proton flow parameter F as a function of the incident energy for Nb+Nb (from Ref. [201]) and Au+Au 



collisions (from Ref. [202]) 
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Fig. 16 shows the average transverse momentum of kaons in the reaction plane as a function of the center-of- 
mass rapidity y cm in Ni+Ni collisions at 1.93 GeV/nucleon and for impact parameters b<4 fm. The dotted curve, 
corresponding to the case without kaon potential, shows that kaons flow in the same direction as nucleons but with 
a smaller flow velocity. The results with only the kaon vector potential are shown by the dashed curve. The kaon 
flow in this case is opposite to that of nucleons, i.e., the appearance of an 'antiflow' of kaons with respect to nucleons. 
With a weak repulsion due to both scalar and vector potentials, it is seen that the kaon flow, shown by the solid 
curve, is in the same direction as that of nucleons but its strength is significantly reduced. Thus, the repulsive kaon 
potential tends to make kaons flow away from nucleons, and its effect depends sensitively on the strength of the kaon 
potential. It is therefore possible to study the kaon potential in nuclear medium by measuring kaon flow in heavy-ion 
collisions. Preliminary data from the FOPI collaboration |57| , shown by the solid circles with error bars, seem to be 
consistent with the existence of both scalar and vector potentials for a kaon in nuclear medium. 
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Fig. 15 Pion average in-plane transverse momentum as a function of rapidity from Au+Au collisions at 1 
GeV/nucleon and impact parameters b=3 (dotted), 6 (dashed), and 9 (solid) fm. (from Ref. |5^| ) 

With the confirmation of potential effects on kaon flow in heavy-ion collisions, it is interesting to note that the 
attractive antikaon and antiproton potentials in nuclear medium are expected to change the antiflow of these particles 
caused by their shadowing from nucleons as a result of large absorption cross sections. Indeed, it has been recently 
shown that including the attractive potential given by the chiral Lagrangian, the strong antikaon antiflow is seen to 



change into a weak flow [298]. Similar medium effects have been shown for pion flow in the QMD calculations with 



an in-medium pion dispersion relation [299] 
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Fig. 16 Kaon average in-plane transverse momentum as a function of rapidity from Ni+Ni collisions at 1.93 
GeV/nucleon. The experimental data from the FOPI collaboration |57| are shown by circles. 

B. heavy-ion collisions at AGS and SPS energies 

Ultra-relativistic heavy-ion collisions have been carried out at AGS with a beam energy of about 10 GeV/nucleon and 
at SPS with an incident energy of about 200 GeV/nucleon. The main motivation for carrying out these experiments 
at increasingly higher energies is to create and study in the laboratory the quark-gluon plasma, which is believed 
to exist in the early evolution of the universe at about one micro second after the Big Bang. At these energies, 
the reaction dynamics becomes much more complex than that at SIS energies, as many hadrons, such as eta, kaon, 
antikaon, rho and omega as well as other higher baryon resonances, are abundantly produced, and we thus need to 
include them explicitly. Since the interactions among these hadrons are not well understood, various assumptions 
have been introduced to relate the unknown ones to the known ones. Nevertheless, cascade-type models have been 
developed and seem to describe the global features of heavy-ion collisions at both AGS and SPS energies reasonably 
. In this subsection, we shall discuss some selected topics such as the strangeness enhancement, the 
difference in the slope parameters of K + and K~ momentum spectra, and the enhancement of very low momentum 
kaons. 



well te7, 181,1* 



36 



1. Kaon enhancement 



In heavy-ion collisions at AGS energies, the K + /tt + ratio has been found to be enhanced relative to that from 
proton- nucleus collisions J63| , |64[ . Different mechanisms have been proposed to explain this enhancement. In RQMD 
[ 182| , p00[ , it is attributed to meson-baryon interactions in the collision. In ARC |18l| ] , it is due to interactions between 
baryon resonances created in the early stage of the collision. In a most recent relativistic transport model (ART) fl97| , 
significant contributions are found not only from the baryon-baryon and meson-baryon interactions but also from the 



meson-meson interaction. On the other hand, it has been shown in Ref. [301] that the measured kaon yield can be 
explained by a simple fireball model at both thermal and chemical equilibrium. Since all models predict a more or 
less equilibrated system at freeze out, it is thus difficult to verify from the kaon yield which models correctly describe 
the kaon production mechanism. Furthermore, in these models the energy density reached in the early stage of the 
collision when most kaons are produced have been found to exceed the one believed to be for a quark-gluon matter. 
The conclusions drawn from these models are thus expected to change when a more proper treatment of the initial 
stage is included. 




Fig. 17 Particle rapidity distributions in central Si+Au collision at 14.6 GeV/c/nucleon. The left, middle, and 



right panels are from Refs. roll], [181], and [300], respectively. The experimental data from Ref. p4J are also shown. 



To understand the enhanced kaon production, the relativistic transport model has been used to study the expansion 
stage of heavy-ion collisions at these energies by assuming that a fireball is formed in the initial stage |^,^l| . The 
fireball then expands, and kaons and antikaons are produced from baryon-baryon, meson-baryon, and meson-meson 
interactions. Both kaons and antikaons with in-medium masses propagate through the hadronic matter under the 
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influence of the mean-field potentials and undergo collisions with both nucleons and pions. Furthermore, antikaons 
can be destroyed via the reaction KN — > nY. 

In the left panel of Fig. 17, the particle rapidity distributions from Ref. |nj are compared with the experimental 
data from Ref. p3[ |. One sees that good agreements with the experimental data are obtained for pions, kaons, and 
antikaons. The failure of the calculated proton distribution at smaller rapidities is due to the neglect of protons from 
the target spectator in the fireball approach. In this model study, a significant number of kaons are produced from 
meson-meson interactions which become important as a result of reduced total kaon-antikaon masses in dense matter 



[136]. 



In the middle panel of Fig. 17, the results from the ARC model of Ref. |181| are compared with the experimental 
data for proton, ir + and K + . Medium effects are neglected in this cascade- type calculation. Kaons are found to 
be mainly (about 70%) produced from the interactions between baryon resonances. The results from the RQMD 
calculation f30(| for 7r + , K + , and K~ are shown in the right panel of Fig. 17. Again, the enhanced kaon production 
can be well described. In the RQMD, kaon production from baryon- meson collisions is found to be important. Similar 
results have been obtained in the ART calculation of Ref. |)7| . The successful description of the kaon data by cascade- 



type calculations, however, does not rule out the mechanism proposed in Ref. |60|. If the mass of KK decreases as 
a result of chiral symmetry restoration, more kaons are produced, but they will be destroyed by reverse processes 
to maintain chemical equilibrium. One thus needs to look at other observables, such as the phi meson yield as it is 
sensitive to the change of KK energy as a result of the small difference, about 30 MeV, between the phi meson mass 
and the sum of K and K masses. The dilepton spectra from kaon-antikaon annihilation will also be useful, as it shows 
directly the invariant mass of the kaon-antikaon pair through the phi meson peak. 



2. Difference in the slope parameters of kaon and antikaon transverse mass spectra 



In Ref. 1 61 1, the transverse mass spectra of particles have also been calculated in the fireball model, and they 
are shown in Fig. 18. All particles have essentially exponential distributions. The slope parameters for kaons and 
antikaons show the difference expected from mean-field effects, i.e. the effective temperature of antikaons is lower than 
that of kaons. This can be understood as follows: Consider antikaons with high initial energies in the fireball that is 
formed in the collision. They move relatively fast and escape thus from the fireball while its size has increased only 
slightly. But these antikaons must use up much of their kinetic energies to climb out of the deep potential well so that 
the measured kinetic energies are substantially smaller. On the other hand, low energy antikaons stay in the fireball 
and escape later in time when the fireball size is large and the mean-field potential becomes negligible small. These 
antikaons therefore do not loose much energy. The net effect is that the apparent temperature of antikaons after freeze 
out is lower than the initial temperature. For kaons, the potential is only slightly repulsive because of cancellation 
between the attractive scalar field and the repulsive vector field. The change of kaon apparent temperature is thus 
very small. If antikaons and kaons have similar initial temperatures, then one expects to see a lower final temperature 
for antikaons than for kaons. 
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Fig. 18 Particle transverse energy spectra from central Si+Au collisions, (from Ref. |pl[) 

Similar effects are seen for antiprotons. The deep attractive mean- field potential for antiprotons makes their appar- 
ent temperature much lower than that for protons ]6^ |. All these results seem to be supported by the experimental 
observations. 

We note that the above discussions are for Si+Au collisions. For Au+Au collisions at AGS energies, the initial 
compression is expected to be much more appreciable than in Si+Au collisions, so a larger radial expansion may 



appear and thus reduces the effects of mean- field potential on particle spectra [286] 



3. cool kaons 

Another interesting preliminary experimental observation in heavy-ion collisions at AGS energies is the enhancement 
of low transverse mass (m t ) K + and K~ in Si+Pb collisions at 14.6 (GeV/c)/nucleon by the E814 collaboration [|302|| . 
The spectra of these extremely cold kaons can be characterized by an inverse slope parameter (temperature) as low 
as 15 MeV, which is about one order-of-magnitude smaller than the temperature of normal kaons measured in earlier 
AGS experiments These low mt kaons cannot be obtained in conventional models such as RQMD Jl82|, ARC 



[181 1, and ART WA . Although this experimental result has not been unambiguously confirmed by other experiments 



031] , it has already generated some interesting theoretical explanations. 
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Fig. 19 Initial (dotted) and final (solid) kaon and antikaon transverse energy spectra. The experimental data are 
from pi. (from Rcf. pi) 



Based on the Georgi vector limit of the hidden gauge theory, Brown and Rho ]304| have argued that the vector 
interaction vanishes when the chiral symmetry is restored. This seems to be supported by calculations based on 
effective hadronic models for the vector coupling constant in medium [)305| . If this is indeed the case, then one 
expects that the kaon potential will become attractive at very high density and temperature as the scalar attraction 
becomes dominant. This idea was first adopted by Koch |306| , who considered the temperature effects, to explain the 



appearance of these low energy kaons. In Ref. [ 307 1 , the density effect is emphasized. Because of the reduced vector 
potential, the nuclear equation of state becomes softened, leading to the possible formation of a density isomer. As 
a result, the initial expansion of the nuclear matter is relatively slow so that kaons and antikaons can be effectively 



cooled by the attractive mean-field potentials. This effect has been studied in the relativistic transport model [307] 
by assuming that the system is initially in the density isomer. The results are shown in Fig. 19 and compared with 
the preliminary E814 data. For both kaons and antikaons, one obtains in the final transverse mass spectra a cold 
low-mt component as observed in the experiment. 



4- strangeness enhancement in heavy-ion collisions at SPS energies 

In heavy-ion collisions at SPS energies, enhanced production of strange particles have also been observed. For 
example, the antilambda yield in the NA35 experiment of S+S at 200 GeV/nucleon [|73[ is 1.5 per event and is 
115 times greater than that in p-p collisions at same energy. Compared with the 36-fold enhancement of negatively 



40 



charged particles, most of them being negative pions, there is a factor of three enhancement of antilambda yield in 
these collisions. This enhancement can be explained if one simply assumes that a quark-gluon plasma is formed in 
the initial stage of the collisions. Other explanations have also been proposed. Aichelin and Werner have emphasized 



the importance of many-body cluster effects [ [308[ . Sorge et al. [309] have shown that the formation of a color rope 
from string excitations can lead to enhanced production of antilambdas. In Ref. j72|, this enhancement is explained 
by the lower antilambda production threshold as a result of reduced antilambda mass in dense matter. In a simplified 
hydrochemical model, the process KM — » AiV, where M denotes either a pion or a rho meson, gives an enhanced 
antilambda production in the collision. A similar explanation based on the relativistic mean field theory has been 
proposed in Ref. ]310[ . 

The reduced phi meson mass in medium has also been shown to lead to an enhanced phi meson production in 
CERN heavy-ion collisions 0,0. 

C. dilepton production in heavy-ion collisions 

Since dileptons are not subject to the strong final-state interactions associated with hadronic observables, they are 
the most promising probe of the properties of hot dense matter formed in the initial stage of high energy heavy- 
ion collisions. Dileptons have thus been proposed as useful observables for studying medium modifications of the 



pion dispersion relation at both finite density and temperature [311,312], the in-medium properties of vector mesons 



p8,313- 316| , and the phase transition from the hadronic matter to the quark-gluon plasma J3 IT - 32 1 1 . 

There are many sources for dilepton production in hadronic matter. These include the proton-neutron and pion- 
nucleon bremsstrahlung, the Dalitz decay of pions, etas, and deltas, the direct decay of vector mesons, as well as 
the pion-pion and kaon-antikaon annihilation that proceed through the formation of vector mesons. To study the 
in-medium properties of vector meson properties, the pion-pion and kaon-antikaon annihilation are most important 
due to vector dominance in the pion and kaon electromagnetic form factors. Neglecting lepton masses, the dilepton 
production cross section from the annihilation of two pseudoscalar meson P through a vector meson V is given by 

8ira 2 k 

3M 6 (M z — my) 1 + (mvtv) 
In the above, M is the invariant mass of the lepton pair, a is the fine structure constant, and k is the three-momentum 
of the pseudoscalar meson in the center-of-mass of the vector meson. The mass and width of the vector meson are 
denoted by my and Ty, respectively. The value of a is 1 and 1/9 for 7r + 7r~ — ► p° — ► e + e~ and K + K~ — > tft — > e + e" 
(K°K° — > (j) — > e + e~), respectively. In this "form factor" approach, the formation time of the dilepton is neglected and 
the effect of the intermediate vector meson is included through the electromagnetic form factor. This is a reasonable 
approach for the case without medium effects. 

With medium-dependent vector meson masses, one needs to adopt the "dynamical approach" by including explicitly 
the formation, propagation and decay of the intermediate vector meson so that the change of its properties in medium 
can be included. In this approach, the vector meson formation cross section from the meson-meson annihilation is 
given by 



6tt (JWT 



"rr^v h ^ {M ^Ziy V +{ m yT v r (66) 
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where & is 2 and 1 for tt + it — > p° and K + K — > <fi (K°K° — > </>), respectively. The decay width of a vector meson 



of mass M to the dilepton is given by [ 142 1 



4na 2 m 



4 



where the coupling constants are gp-n^l^n « 2.9 and g^Kx/^ ~ 1-7, determined by their respective decay widths. 

We note that in the dynamical approach the vector meson mass and width in Eqs. ( |66| ) and ( |67| ) are evaluated 
using the local density at the time of its formation and decay, which may be different due to the change in density. 
It has been found in Ref. |p8| that for vector mesons with large decay widths (e.g., the rho meson), the form factor 
approach gives results qualitatively similar to the dynamical approach. However, for mesons with small decay widths 
(e.g., the phi meson), the dynamic approach is required as the medium effects are overestimated in the from factor 
approach. 




Fig. 20 Dilepton invariant mass spectrum from Ca+Ca collisions at 1.05 GeV/nucleon. The experimental data 



from the DLS collaboration [322] are given by solid squares while the theoretical total yield is given by the thick solid 



curve, (from Ref. [ 218 1 ) 



In the following, dilepton production from heavy- ion collisions at SIS/GSI, SPS/CERN, and RHIC/BNL is studied 
in the relativistic transport model, and the results are compared with available experimental data. 
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1. dileptons from SIS /G SI 



Dileptons have been measured from heavy- ion collisions at Bevalac by the DLS collaboration [322]. The experi- 
mental dilepton invariant mass spectrum from Ca+Ca collisions at 1.05 GeV/nucleon is shown in Fig. 20. Based 
on transport models, perturbative calculations of dilepton production, including contributions from proton-neutron 
bremsstrahlung, pion-nucleon interactions, delta and eta decays, and pion-pion annihilation, have been carried out in 



Refs. 21S.32S.324]. A typical result from Ref. [218] is shown in Fig. 20. It is seen that dileptons with small invariant 
masses are mainly from proton-neutron bremsstrahlung and delta decays while dileptons with large invariant masses 



arc dominated by pion-pion annihilation. The contribution from eta decays, that was neglected in Rcf. [218], turns 



out to be very important for invariant masses below about 500 MeV [324]. Nevertheless, the eta has a width of only 
1.2 keV and decays thus outside of the matter. Its contribution to dilepton production can in principle be subtracted 
out if its spectrum is also measured. We also see clearly the rho meson contribution to dilepton production from 
pion-pion annihilation as a result of rho meson dominance in the pion electromagnetic from factor. The data are, 
however, not accurate enough to extract the in-medium rho meson properties. 



Ni+Ni, 2.0 GeV/nucleon, b=0 fm " 
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Fig. 21 Dilepton invariant mass spectra from pion-pion and kaon-antikaon annihilation in Ni+Ni collision at 2 
GeV/nucleon and impact parameter b=0 fm. The dashed and solid curves are obtained using free and in-medium 
meson masses, respectively, (from Ref. (58)) 

Future experiments at SIS with the HADES detector |39) will provide data with vastly improved statistics to 
allow us to study more clearly the medium effects on vector meson properties. To see the medium effects in these 
experiments, the dilepton invariant mass spectra from Ni+Ni collisions at an incident energy of 2 GeV/nucleon and 
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impact parameter of fm with free vector meson masses (dashed curve) and in-medium vector meson masses given 



by Eqs. (43) and (|44| ) from the QCD sum rules (solid curve) have been calculated |58j and are shown in Fig. 21. 
For pion-pion annihilation, medium effects shift the rho meson peak to around 550 MeV, and its height increases by 
about a factor of four. For kaon-antikaon annihilation, the phi peak shifts slightly to a lower invariant mass, and its 
width is also broadened when medium effects are included. In addition, there appears a shoulder around 950 MeV. 
For a complete picture of dilepton production at SIS energies, contributions from other channels, such as the Dalitz 
decay, bremsstrahlung, vector mesons produced in initial nucleon-nucleon collisions, and vector mesons from decay of 



high resonances 1 325 ] , need to be included. 

For heavy-ion collisions at AGS energies, a similar shift of the rho meson strength to lower dilepton invariant masses 



has been seen in RQMD calculations using also the QCD sum-rule results for in- medium vector meson masses [326|. 
Since a omega meson has a much longer lifetime and decays mostly at freeze-out, an interesting phenomenon of p-cu 
splitting was seen in the dilepton spectra. It would be very useful if dilepton production from heavy-ion collisions at 
AGS energies could be measured. 



2. dileptons from SPS/CERN 

For heavy-ion collisions at SPS energies, hot and dense matter is also formed in the initial stage of the collisions. One 
expects that medium effects will lead to a shift of the vector meson peaks in dilepton invariant mass spectra. Experi- 
ments from both the CERES (67| and the HELIOS-3 |38| collaboration have shown that there is an excess of dileptons 
over those known and expected sources which cannot be explained by uncertainties and errors of the normalization 
procedures |6^]. In particular, in the CERES experiment for central S+Au collisions at 200 GeV/nucleon, a signifi- 
cant enhancement of dileptons with invariant masses between 250 MeV to 1 GeV over that from the proton-nucleus 
collision has been found. 

These experimental data have generated a great deal of interest in the heavy-ion community. Different models, 



ranging from schematic estimates based on a possible enhancement of i] and rf production [327, 328 1, to hydrodynam 



ical models [329,3301, and relativistic transport models [^,^6,331 333), have been used to study this phenomenon 



Although the contribution from pion-pion annihilation, that has not been included in the 'cocktail' of the CERES 
collaboration, was found to be important for low-mass dileptons, the data in the low mass region are still above the 
theoretical results from these calculations when vector meson properties in free space are used. This situation is 
summarized in ig. 22 where the results from Ref. f66f| (solid curve) using the relativistic transport model with initial 



conditions determined by the RQMD, Ref. [ 331 1 (dashed curve) from the Hadron-String Dynamics, and Ref. [32£] 
(dotted curve) from the hydrodynamical model are shown together with the CERES data. It is seen that in all three 
studies the calculated low mass dileptons are by about a factor of 2-4 below the data, and for dileptons around m pM 
they are slightly above the data. It is also of interest to note that the three calculations, although with different 
dynamical models, agree more or less with each other, and the differences are less than 20%. A strong peak around 



m<£ in the results of Ref. [ 32E ] will become a bump once the mass resolution of the CERES collaboration is properly 
included. 



44 



central S+Au, 200 GeV/nucleon _ 




0.2 0.4 0.6 0.8 1 1.2 
M (MeV) 



Fig. 22 Dilepton invariant mass spectra in central S+Au collisions at 200 GeV/nucleon with free meson masses. 



The solid, dashed, and dotted curves are from Refs. p6j, [331], and 32£], respectively. The experimental data from 



the CERES collaboration |37j are shown by solid circles, with bars denoting the statistical errors while the systematic 
uncertainties are given by brackets. 

The results using in-medium vector meson masses are shown in Fig. 23 by the solid curve. The agreement with the 
experimental data is greatly improved. The enhancement at low invariant masses is due to both the shift of primary 
rho mesons to lower masses and pion-pion annihilation occurring in hot dense matter. Since pions have a thermal 
distribution, most pion pairs are of low invariant mass. When the rho- meson mass is reduced, its formation probability 
from pion-pion annihilation is enhanced, thus increasing the production of low-mass dileptons. The remaining peak 
around to p . w then comes from the decay of omega mesons which have a very small decay width, and therefore mostly 
decay at freeze out when their mass has returned to the free value. One notes that dileptons with invariant masses 
below about 300 MeV are mainly from the Dalitz decay of 7r°, u>, and r\. 

Dimuon invariant mass spectra have been measured by the HELIOS-3 collaboration in S+W collisions at 200 
GeV/nucleon ^38|. The data also show an enhancement of dileptons around M ~ 0.4 — 0.6 GeV. This provides 
another possible indication that the vector (rho) meson mass might be reduced in hot dense medium. The comparison 
of theoretical results obtained from the relativistic transport model with the HELIOS-3 data is shown in Fig. 24 ||S6| . 
With free vector meson masses, the theoretical results are below the HELIOS data in the mass region from 0.35 to 0.6 
GeV by about a factor of two, and slightly above the data around m pM . Using in-medium vector masses, we again 
see an enhanced dilepton production in the low mass region and a reduction around m p ^ . This brings the theoretical 
results in better agreement with the data. Similar results have been obtained in Ref. [)332| using the Hadron-String 
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Dynamics. 
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Fig. 23 Same as Fig. 22. The dotted and solid curves are obtained using free and in-medium meson masses, 



respectively, (from Ref. |66|) 



The medium effects seen in the HELIOS-3 case are smaller than that in the CERES case. This is due to the fact that 
in the HELIOS-3 experiment dileptons are measured at forward rapidity (3.7 < rj < 5.5) where the charged-particle 
multiplicity is low, while in the CERES experiment they are measured in the midrapidity (2.1 < rj < 2.65) with a 
higher charged-particle multiplicity. It is also of interest to note that the theoretical results from both Ref. p6f and 



Ref. [332 1 are well below the HELIOS-3 data above 1.2 GeV. This may indicate the importance of processes such as 
ircii — > l + l~ [334,335] and the decay of heavier vector mesons such as w(1390) [336,337]. Also, in this mass region, 



the contribution from the quark- gluon plasma and the initial Drell-Yan processes might become important. 



3. dileptons from RHIC/BNL 



In heavy-ion collisions at the Relativistic Heavy-Ion Collider (RHIC) being constructed at the Brookhaven National 
Laboratory, dileptons will be measured as they are likely to carry the signature for the quark-gluon plasma expected 



to be created in the collisions [ 317 319 , 321 ] . However, there are a number of different sources of dileptons in heavy- 
ion collisions at these energies. For dileptons with low invariant mass, i.e., below about 1 GeV, they are dominated 
by hadronic processes, such as 7r° and r/ Dalitz decays as well as vector meson (p, u>, and 0) direct decays. For 



intermediate-mass dileptons, i.e., between 1 and 2 GeV, dilepton production from heavy meson resonances [ 335,336) 
and charmed meson (D and D*) decays become important. It is in this invariant mass regions that it may be possible 
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to see dileptons emitted from the quark-gluon plasma via quark-antiquark annihilation if its temperature is sufficient 



high 1 33?]. Above 2 GeV, dileptons are from Drell-Yan process involving the annihilation of quarks and antiquarks 
in the initial colliding nuclei. In this region of dilepton invariant mass, an interesting proposal to see the signature of 



the quark-gluon plasma is the J/vP suppression 339 due to its dissociation in the quark-gluon plasma as a result of 
color screening. Current experiments at CERN/SPS have indeed shown that the number of J/vP observed versus the 



Drell-Yan background is reduced by about a factor of 2 in heavy-ion collisions than in proton- nucleus reactions [340 



Unfortunately, the data can also be explained by hadronic absorption models [341-343] 
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Fig. 24 Dimuon invariant-mass spectrum from pion-pion annihilation in central S+W collisions at 200 GeV/nucleon. 
The dotted and solid curves are obtained with free and in-medium masses, respectively. The experimental data from 
the HELIOS-3 collaboration are shown by solid circles, (from Ref. fl66j ) 

On the other hand, the change of phi meson mass in hot matter as discussed in Sect. II. C has led to a new 
signature for identifying the quark-gluon plasma to hadronic matter phase transition in ultrarelativistic heavy-ion 



collisions [344]. In a boost invariant hydrodynamical calculation with transverse flow and including temperature- 
dependent vector meson masses, it has been found that a low mass phi peak at ~ 880 MeV besides the normal one 
at 1,020 MeV appears in the dilepton spectrum if a first-order phase transition or a slow cross-over between the 
quark-gluon plasma and the hadronic matter occurs in the collisions. Assuming an initial temperature T,- t = 250 MeV, 
a critical temperature T c — 180 MeV, and a freeze out temperature Tf = 120 MeV, the results are shown by the solid 
curve in Fig. 25. The three peaks correspond to the omega, the low-mass phi, and the normal phi meson, respectively. 



Because of their reduced mass at finite temperature according to the QCD sum- rule studies |345 ], decays from rho 
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mesons lead to a broad distribution of dileptons at low masses. 

The low-mass phi peak is due to the nonnegligible duration time for the system to stay near the transition tem- 
perature (about 10 fm/c) compared with the lifetime of a phi meson in the vacuum (~ 45 fm), so the contribution 
to dileptons from phi meson decays in the mixed phase becomes comparable to that from phi meson decays at freeze 
out. Also, its width from collisions with other hadrons remains small (about 10 MeV) §4(||47|]. 

Without the formation of the quark-gluon plasma, the low mass phi peak is reduced to a shoulder in the dilepton 
spectrum as shown by the dotted curve in Fig. 23, which is obtained by assuming that the initial temperature of the 
hadronic matter is just below T c . 




Fig. 25 Dilepton invariant mass spectra from ultrarelativistic heavy-ion collisions. Solid and dotted curves corre- 



spond to scenarios with and without the formation of a quark-gluon plasma, respectively, (from Ref. [ 344 1 ) 



Since the transverse expansion velocity during the mixed phase is relatively small, the transverse momentum 
distribution of the low mass phi mesons is largely determined by the temperature of the mixed matter. The low mass 
phi mesons thus also provide information about the temperature at which the quark-gluon plasma to hadronic matter 
transition occurs. 



V. SUMMARY AND OUTLOOK 

Because of the partial restoration of chiral symmetry, the quark condensate decreases in hot dense matter. According 
to the prediction of QCD sum rules, hadron masses are then reduced in medium. As a result, the thresholds for particle 
production in medium are reduced, and the cross sections for their production are enhanced due to the increase of 
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phase space in the final states. The dense matter also gives rise to strong vector mean-field potentials. With these 
medium effects consistently included in the relativistic transport model, it has been shown that for heavy-ion collisions 
at SIS energies they can account for the observed enhancement of kaon, antikaon, and antiproton yields, the enhanced 
soft pions in the transverse direction, and particle flows in these collisions. 

To place these conclusions on a firmer ground, a number of issues need to be addressed in future investigations. In 
heavy-ion collisions at SIS energies, the colliding system consists mainly of nucleons, baryon resonances, and pions. 
Particles can thus be produced from baryon-baryon, baryon-meson, and meson-meson collisions. Because of the small 
center-of-mass energy of the pion and the fact that pions materialize in the expansion stage of the collisions, hadron 
production from pion-pion interactions can be safely neglected at SIS energies. Kaon and antikaon production from 
pion-baryon collisions have been calculated and found to account for about 20% of their total yield. With improved 
experimental data, these contributions need to be more carefully studied. Similarly, antiproton production from 
pion-baryon collisions has not been considered and certainly need to be evaluated. 

For subthreshold particle production, the most important contribution comes from baryon-baryon interactions, of 
which interactions involving baryon resonances are particularly important. One needs thus particle production cross 
sections in baryon-baryon collisions. For proton-proton collisions there exist empirical data for kaon, antikaon, and 



antiproton production, mostly at high beam energies 1 348 ] . The general strategy in the transport model is to introduce 
a parameterization for these data and then to extrapolate it to the threshold. The production cross section near the 
threshold then depends on the way the parameterization and extrapolation are introduced. This is particularly serious 
for antiprotons since there are essentially no experimental data below 20 GeV. Moreover, no experimental data exist 
for particle production from baryon-baryon collisions involving baryon resonances. Various assumptions have been 
made to relate these cross sections to that in proton-proton collisions. The validity of these assumptions certainly 
need to be carefully assessed. Also, the role of baryon resonances that are higher than the delta needs to be further 



investigated |49j , g50 |. 

For antikaon and antiproton, there is an additional complication from their large annihilation cross sections. This 
raises the question of numerically how to accurately treat the annihilation process in transport models. Indeed, 
different treatment of antiproton annihilation could lead to almost an order of magnitude difference in the final 



antiproton survival probability [51,52 241 242,351]. Moreover, the antikaon and antiproton annihilation cross sections 



might be modified in medium due to changes in their masses, the screening effects from mesons [181|, and the effect 



of Bosc-Einstein enhancement in the final state [352]. For antiproton annihilation, which leads to about five pions, 
the latter effect is expected to be important in heavy-ion collisions at AGS and SPS energies as the abundance of 
pions is large in these collisions. 

The medium modification of hadron masses has so far been included in transport models only in a minimum way 
by changing the threshold of particle production. In terms of Feynman diagrams, this amounts to change only the 
masses of the initial and final hadrons. In principle, the properties of exchanged particles can also be modified p0,231|. 
Besides changes of hadron masses, the coupling constants and the form factors in the vertices of Feynman diagrams 
might be modified in medium as well. All these need to be consistently studied. 

Finally, off-shell propagation of particles has been neglected in all transport models. This effect on subthreshold 



particle production has been shown to be small [353]. However, off-shell effects due to correlations in the initial wave 
function may be important and need to be further studied. 
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For heavy-ion collisions at AGS energies, medium effects can also explain the enhanced K + /n + ratio, the difference 
between the slope parameters of K + and K~ transverse kinetic energy spectra, the lower apparent temperature 
of antiprotons than that of protons, and the cool kaons with low transverse masses. However, these studies have 
been based on a fireball initial conditions so the initial noncquilibrium dynamics is not included. There are already 



models, such as the RQMD [ 182 1 , ARC |181|, and ART [|97[ , which can treat the initial stage of heavy-ion collisions at 
these energies as enough inelasticity has been introduced in these models. It will be useful to extend the relativistic 
transport model by including these additional processes so both the collision dynamics and medium effects can be 
properly included. 

Medium effects are also relevant for heavy-ion collisions at higher energies from the SPS at CERN. The reduced 
particle production thresholds in medium have been found to lead to enhanced production of phi mesons, antilambdas, 
and low-mass dileptons. As in the case of heavy-ion collisions at AGS energies, the relativistic transport model at 
its present form can only describe the expansion stage of these collisions. In the initial nonequilibrium stage, particle 



production from string fragmentation has shown to be essential [182]. It is thus again important to extend the 
relativistic transport model to include consistently in the initial stage the string dynamics. Such a step has already 



being taken by Cassing et al. [ 354 1 . 

We have shown in this review that high energy heavy-ion experiments offer the possibility to study the properties 
of hadrons in the hot dense matter formed in the initial stage of the collisions. This study is not only of interest in 
its own right but also important for future heavy-ion experiments at ultrarelativistic energies where the quark-gluon 
plasma is expected to be formed. To find the signatures of the quark-gluon plasma, it is essential to have a good 
understanding of the hadronic matter that exists both in the initial and final stages of heavy-ion collisions. Moreover, 
the modification of hadron properties in hot dense matter is expected to give rise to new signals for the quark-gluon 
to hadronic matter transition. For example, a low mass phi peak besides the normal one has been predicted to 
appear in the dilepton invariant mass spectra from ultrarelativistic heavy-ion collisions as a result of the decay of phi 
mesons with reduced in-medium mass during the transition. Furthermore, measurements of the transverse momentum 
distribution of these low mass phi mesons offer a viable means for determining the temperature of the transition. 

However, experience from heavy-ion collisions at SIS, AGS, and CERN energies has shown that to extract useful 
physics from heavy-ion collisions requires transport models which can describe properly the entire collision dynamics. 
The same is expected for heavy- ion collisions at RHIC, so it is important to have a model which includes the initial 



parton cascade |35q| , the evolution of the resulting quark-gluon matter [356], the dynamics of hadronization [ 357 1 , 
and the final hadronic dynamics. In this respect, the relativistic transport model developed extensively for low energy 
heavy-ion collisions provides a link of future developments in ultra-relativistic transport model to the experimental 
observations. 
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